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di ; dil 


veometrisches Wegelement ; n = 


chungszahl (des « Mediums ») an der Stelle von dl]. 
Dementsprechend bezieht sich auch das Prinzip der 
Konstanz der optischen Weglingen > (nr. Al) bzw. 


2P. 

| n dl fiir zwei oder mehr Strahlen zwischen zwei 
Pp 

e 1 


Punkten P, und P, zunachst nur auf benachbarte 
Strahlen. Bei seiner Anwendung auf « aplanatische » 
Flichen oder auf « cartesische » Flachen — _ aiso 
Flachen, bei denen nur ¥ (n. A 1) = const fiir alle von 
einem Punkte ausgehenden Strahlen zwischen diesem 
Punkte und seinem Bildpunkte gefordert wird, wah- 
rend ja bei den aplanatischen Flachen ausserdem fiir 
alle diese Strahlen die Erfiillung der Sinusbedingung 
erforderlich ist —- handelt es sich indessen um eine 
Vielzahl, eine Mannigfaltigkeit von Strahlen, die ‘ei- 
neswegs alle als « einander benachbart » zu bezeichnen 
sind, wohl aber besitzt hier jeder der betrachteten 
Strahlen Nachbarstrahlen, mit denen er bzw. seine 
optische Weglainge verglichen und als optisch 
gleich lang verlangt — wird derart, dass diese « Gleich- 
heit » iiber das ganze Strahlenbiindel fortgesetzt wer- 
den kann, ohne die Forderung der « Nachbarschaft » 
bei den einzelnen Vergleichen aufzugeben. 


|. Unterschied der optischen Weglangen eines von 
einem Achsenpunkt ausgehenden, zur Achse (beliebig) 
geneigten Strahls und des zugehorigen axialen Strahls 
zwischen dem gemeinsamen Objektpunkt und ihrem 
bildseitigen Schnittpunkt bei beliebig grosser sphari- 
scher Aberration. — Bei der Berechnung der optischen 
(Licht-) Wege zweier beliebiger, vom gleichen Objekt- 
punkt ausgehender Strahlen von diesem Objektpunkt 
bis zu ihrem bildseitigen Schnittpunkt bzw. bei der 
Berechnung ihres optischen Wegunterschiedes handelt 
es sich aber um diskret liegende Strahlen — z. B. 
een oder mehrere (bzw. unendlich viele) Strahlen 
bestimmter Neigung zu dem « Hauptstrahl » des abbil- 
denden Strahlenbiindels und diesen « Hauptstrahl » 
selbst —, ohne dass hier die Berechtigung der Anwen- 
dung des Fermatschen Prinzips oder eines ihm ana- 
logen besteht, da es sich jetzt ja nicht — wie bei den 
aplanatischen oder cartesischen Flachen — um eine 
kontinuierliche Mannigfaltigkeit von « benachbarten » 
Strahlen handelt. 


Es sei nun zunachst fiir den Fall der spharischen 
Aberration der Unterschied der optischen Weglangen 
eines zur Achse geneigten und eines mit der Achse 
zusammenfallenden; vom gleichen Objekpunkt ausge- 
henden Strahls in ihrem bildseitigen Schnittpunkt 
berechnet. 


Dabei gehen wir aus von einer fritheren Arbeit des 
Verfassers [1], in der die Gleichungen fiir die Wellen- 
flachen eines Strahlenbiindels abgeleitet wurden, die 
zu einem Strahlenbiindel mit spharischer Aberration 
gehoren. 


Die spharische Aberration, die wir wie iiblich durch 
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\s’ bezeichnen, lisst sich (in allen Fallen) darstellen 
durch die (im allgemeinen endliche) Reihe 

(1,1) As’ =a tg? u’ + btgtu' + c tg* w’+ dtg*u’ +..., 
wo u’ der — wie itblich gemessene — Winkel der 
(bildseitigen) Strahlen des Biindels gegen seine Sym- 
metrieachse (Achse des optischen System) sei. 

Wie in der angegebenen Arbeit gezeigt wurde, lau- 
ten die Gleichungen der zugehérigen Wellenflachen 
dann — bezogen auf ein Koordinatensystem, dessen 
2-Achse mit der optischen Achse zusammenfallt, 
(1,2) €=—2a'+a’ tg?u’ +0’ teu’ + ec’ te we 

+... — const. cos w! 


2 4 6 
(13) %pe= if a’tgu'’ + x bite 5 An eA Teer 


== COnsta Silmraas 


wenn e den Abstand der Punkte der Wellenflachen 
von der Achse bedeutet. Hierin ergeben sich die 


a’, b’, c', ... aus a, b,c, ... des As’ durch die Formeln 
vp 
—— mae paee VV 
a a 3 
6 
b' = b— see G 
(1,4) 
ea ee manny |: 
C c - 
10 
d! = d as ey e! 
lassen sich also fiir praktische Falle — in denen die 


Reihe fiir \s’ abbrechen wird — stets leicht berech- 
nen. In Abb. 1 und 2 sind fiir zwei Strahlenbiindel mit 
spharischer Aberration diese Aberration, ferner Strah- 
len, Wellenflachen und rotationssymmetrische sowie 
entartete Kaustikschale dargestellt. 

Fiir unsere Zwecke des Vergleichs der Linge der 
Lichtwege der die Achse in einem (beliebigen) Punkte 
treffenden Strahlen mit dem Lichtweg des Achsen- 
strahls selbst bis zu Jenem Punkte benutzen wir keine 
Voraussetzung liber das Abbrechen oder Nichtabbre- 
chen der Darstellung von As’ (wohl aber natiirlich bei 
einem Nichtabbrechen die Voraussetzung der Kon- 
vergenz der Reihe). 

Tragen wir vom Achsenschnittpunkt eines der ausser- 


axialen Strahlen — etwa des bildseitigen Neigungs- 
winkels uw, gegen die Achse — eine Strecke (bei 
Asi > 0 : entgegen der Ausbreitungsrichtung des 


Lichtes) ab, die gleich dem Betrage As, der sphari- 
schen Aberration des betreffenden Strahles ist, so 
gilt fiir die den betreffenden Strahl dort orthogonal 
schneidende Wellenflache nach (1,3) 


1 
7; Za’ 4 
COs U, 


y 
. / 
sin W, 


(1,5) = As, = 


4 / 
b tg? = 
3 gr uy 


6 ' 
ae ee tg* u, + ...) + const, 


PT a ee re ee ee Oe ee 
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so dass die zu dieser Wellenflaiche gehérige Konstante 
den Wert hat 


(1,6) 1 4 


(2a’ + 


I 
const, = As} b! tgzu, + 


tf 
cos U, 


6 ' 
ai =a ELSI ree 
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&-und oe -Werte der einzelnen Punkte der betreffenden 
speziellen Wellenfliche liefert, die vom Achsen- 
schnittpunkt des Strahls w, lings dieses Strahls den 
Abstand — As, besitzt. 


Fir uw’ = 0 wird (9;)) = 0 — wie erforderlich — und 


I 4 te? wi 
Mit dieser Konstanten lauten die Gleichungen der (Gre — 2 a ee ae ees . 
betreffenden Wellenfliche Us Rte EN 
- ae COs uw’ \ 6 tg*u, ; 
ao. 2 yt «DoE a eet 7 See NK, 
\ COS U, / 5 cos u, 
/ Pe Se Setzen wir hier fiir \s, den Wert aus (1,5) ein, so wird 
It 2 tf y I L i y 
Bee 08. Wie tb COB L ie okie) 
\ COs U, 
/ 6 foal (1,7) (E10 = — 2a’ — const, . 
| ° i, . t 
= 2 Fh CEG oe zs ~ COS w' | + Da andererseits fiir den Strahl mit vu, = 0 auch 
: S : I : ; 
Oe uet \s, = 0 ist und die entsprechend wie oben kon- 
+ +. = ds, . COS UW’ — struierte Wellenflaiche durch die Kaustikspitze geht, 
tou! 
O5 os 
a 
oO 
a4 oe 
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gf 
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Ass, 1. — Strahlenaberration, Strahlenverlauf, Wellen flachen und Kaustik eines mit spharischer Aberration 
behafteten Strahlenbtindels fiir den Fall, dass 
As’ = atg? u’ mit a = 1000 (baw. =) 
) 
Sys if Ge sin u’ \ so gilt fiir diese Wellenflache nach (1,6) consts = — 2 a’ 
SOE Oh Ra ea aan Sie und (&)) = 0, wobei also der Koordinaten-Ursprungs- 
\ ; ' punkt in die Spitze der Kaustik gelegt ist. 
hats UBS ean oc \ Sollte also die Wellenflache, die den Strahl der 
= Bl ees S —= 4 (ay 5 
a a b' | tet ras sin Achsenneigung tg u, im Abstand — As, von seinem 
, : Achsenschnittpunkt trifft, durch die Spitze der Kaustik 
ieee ieee tg*u, . healed gehen, sollten also — anders ausgedriickt — Achsen- 
-—C ie) — sin u {bse Sn OX Be x Se NTE eee tae ee: 
5 eas w ] strahl und Strahl der Neigung tg wu, in ihrem Schnitt- 


worin uw! wieder der variable Parameter ist, der die 


punkt gleiche Lichtweglingen vom genreinsamen, auf 
der Achse liegenden Objektpunk aus zuriickgelegt 


tis 
4 


haben. so miisste auch const; = — 2 a’ sein, was 
indessen nach (1,6) nicht der Fall ist. Der Unterschied 
der beiden Konstanten ist demnach der Lichtwegun- 
terschied, den der achsengeneigte Strahl gegen den 
Achsenstrahl in ihrem gemeinsamen Schnittpunkt 
besitzen. Bezeichnen wir diesen Wegunterschied durch 


/ 
A Wij9 = Ww (tg u,) — Wy , 


so ist A Wy = (Es)o — (E)o = (1) , also 
ie L te / 
AW) =—2¢—aAs, + ——7 | 2a'+>..0' te? a, 
COS U, 3 
+ —c' tgtu,+ 


0 As! 
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. if os 
bei spharischer Aberration, nur gilt, wenn uv, geniigend 
klein ist, so dass man beide Strahlen als (unendlich) 
benachbart bezeichnen kann. Denn (nur) fiir diesen 
Fall wird Ja A 4). = 0. 


I1. —- Gleichungen des Meridianschnitts der Wellen- 
flachen eines komabehafteten Strahlenbiindels und 
Unterschied der optischen Weglangen eines zum 


Hauptstrahl geneigien Strahls der Meridianebene und 
des Hauptstrahls zwischen dem gemeinsamen Objekt- 
punkt und ihrem bildseitigen Schnitipunkt. — [nt- 
sprechend zu der Darstellung der spharischen Aberra- 
tion in ihrer Abhiingigkeit vom Tangens des bildsei- 
tigen Neigungswinkels w’ der einzelnen Strahlen des 


50 


95\e ds‘ 100017'u'— 4000 1g"u' 
een Vemed X\ , 


a b 
\BB. 2. — Strahlenaberration, Strahlenverlauf, Wellenflachen und Kaustik eines mit spharischer Aberration 
behafteten Strahlenbiindels fiir den Fall, dass 
1 2» 
As’ = atg2u'’ + b tg4u’ mit a = 1000 (baw, ce ): b = — 4000 (va. — = 

) o 
1 \ / A tg? u, abbildenden Strahlenbiindels zur Symmetrieachse des 
SS eed) (2 a’ + A s, | — ——— xX abbildenden Systems lisst sich auch’ der Komafehler 
COR, 7 =x COS, als Potenzreihe der Neigungswinkel x’ der Strahlen 


I i 
x (a’ + 0 tg? un, + ¢’ tg*u, 4..). 


Aus diesem Ausdruck ersieht man also auch unmit- 
telbar, dass der Satz von der konstanten Lichtweg- 
linge zweier vom gleichen Objektpunkt ausgehender 
Strahlen von diesem Objektpunkt bis zu ihrem gemein- 
samen Schnittpunkt nicht in dieser Allgemeinheit 
gilt, auch nicht, wenn der eine der beiden Strahlen der 
Achsenstrahl ist, sondern dass dieser Satz auch hier, 


des Strahlenbiindels gegen dessen Schwerstrahl (Biin- 
delachse) darstellen, sofern wir uns auf Strahlen des 
Meridianschnittes beschrinken. 


Bezeichnen wir — als Analogon zu As’ — den Ab- 
stand des Schnittpunktes eines der im Meridianschnitt 
verlaufenden Strahlen mit dem Schwerstrahl vom 
Schnittpunkt des Schwerstrahls mit einem ihm unmit- 
telbar benachbart verlaufenden Strahl des Meridians- 
chnittes durch Al, so lasst sich Al stets als endliche 
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oder unendliche Potenzreihe von tg x’ darstellen : 


Cae L WEA Se gh tea =P a, tee Aas bee sche 
BGs Oar tno. 
oder, mit tex’ = es 
Y 
a dl a d, FAs 
2,2 A eaetee ee Ol 
(2,2) ins ae aa asl 


Die von Mertrsé (W. Merré, ZS. f. Phys. 33, 1925, 
933) zur Berechnung der Wellenflachen bei sphirischer 
Aberration des betreffenden Strahlenbiindels angege- 
benen Formeln 


. y y o (y) 
V1 a= y e/ const V1 oe Of 
5 4 y 0 
Ca : ( 12) dye 
V1 == y? e/ const V4 + OP 
er 
Sl eee eee 
| Yee Sa ie 
mit 
{late Mater ite Be { HE 
Dy 4 co) = ct = 
ee et ae) 
behalten auch fiir den Komafehler —- wenigstens fiir 


die uns zunichst interessierenden Verhaltnisse im 
Meridianschnitt — ihre Giiltigkeit, wenn wir jetzt die 
2-Achse mit dem Hauptstrahl des Biindels zusammen- 
fallen lassen und die e -Achse senkrecht zur € -Achse 
in der Meridianebene liegt. Man iiberzeugt sich leicht 
davon, wenn man aus (2,3) die Tangente an die Kurve 
(€, e) in irgend einem Punkte (y) der Kurve berechnet. 
Man erkennt dann, dass diese Tangente senkrecht zu 
dem Strahl des (ebenen) Strahlenbiindels legt, dessen 
Neigung gegen den Schwerstrahl (Hauptstrahl) des 


4 
Biindels durch tg x’ = Fa bestimmt ist. 


Zur Berechnung der Wellenflachen — oder besser : 
der Wellenlinien (in der Meridianebene) des komabe- 
hafteten Strahlenbiindels bilden wir nach (2,4) zu- 
nachst 


ay he = Choy ls 
OD (y) = —— 2 oes 
y ] 7] 
(ig 
RO ad eee = 
Os 


*y D(y) dy 
dy = — 2m | = 
4 US 


e’ const V1 -|- We 


LD > 
XS. dy 
~ ed ya 1 Ndoy— 1 | EAA raat 
Uf =r \/ | + y’ 


22) i dy 
oy (eet A) a, | = rn rah 
‘ Y 
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Wir bezeichnen die in der ersten S auftretenden 
Integrale durch /,,,4,, die in der zweiten » auftre- 
tenden Integrale durch /,),,,;.Bei den J.,,4,) handelt 
es sich um die gleichen Integrale, die auch bei der 
Berechnung der Wellenflichen eines mit spharischer 
Aberration behafteten Strahlenbiindels auftreten. Fiir 
sie gilt 

{/n 


n \ 
Ny 


Vi | — 7S = 
(2+ Wefan aS (—1)4 a ao 
cs (2 aoa 

Ty) 


KEntsprechend findet man durch wiederholte par- 
tielle Integration und dadurch gewonnene Zuriick- 
fihrung der Integrale auf niedriger indizierte 


ape Sef! ye 


ap Dont = 


= (—1)nr el Ar Sin a ‘ 


Jen \ n / y 
/9) ‘ = 
ae a ] 5) (—1)4— 7x ~ = < 
3 / 2(n— ? 
u=1 AU | Plats ) (20-2) +1) 
vity 
yrt+1—s) 
Wir erhalten so * 
‘4 ®O 
| UE a 
iN, aiagy® 
a ( in A 
“ Ht | 4/2 \ ae \ ) 
head dan ——# > ( 1) - ‘aa 
n=1 |u=o ane 
\ ee 
~ v1 sale / 2n 
a x Can—1i 4 a(n-+1) n\ n } 
\ N ! yr aie 


eh = , 
ua gu ( 2") ) (on — wy 4-1) 


\ a—yp 
n An 
Qen—1\ n ) 3 


+ { Xr Sin — | SS 


Xs / 


(— bye (2n—1 


hb 


n=1 


Schreiben wir die hier auftretenden Summen aus- 
fuhrlich, so lassen sich wieder — wie bei der Behand- 
lung der spharischen Aberration — gewisse Zusam- 
menfassungen durch Einfiihrung neuer Koeffizienten 
durchfiihren. 

Wir erhalten 

(2,6) | eee dy = 

JVi+y 


n 
(*) Statt nt 


U 
a ate (Bes 


‘ecria kee, 


lasst sich auch schreiben : 


6 Js 


{ + 72 B, 
Very | Oo Bah ek ee ee ee 
y y y" y® y® 


B, Wt Sin = 

y 

Die hier auftretenden Koeffizienten B;(mit « = 1, 2,...) 
ergeben sich aus den Koeffizienten der Reihenent- 
wicklung (2,1) oder (2,2) von A/’ nach folgenden 
Rekursionsformeln : 


= 2. ee 2 : 
Bre, 3 By | 5 
6 | 5 
Bi=a—zB, | B= a— Be 
8 ie 
Bae de — Bs Bs == th, 5 Be 
7) / ) 
10 ) 
Be=ae—— By | B,=a,— 5B, 


Bricht die Reihenentwicklung (2,1) oder (2,2) von 
Ml’ ab, so brechen auch die Koeffizientenreihen der 
Bj automatisch ab. 

Aus (2,6) und (2,3) ergeben sich so die Koordinaten 
der Schnittkurven der Meridianebene mit den Wellen- 
flachen des komabehafteten Strahlenbiindels : 


& = —2B, + (B, tg? x’ + B,tg* x’ + 
\ +- B, tg® x’ + B, tg*® x’ +... 
(2,8) + (B, tg x’ + B, tg® x’ + B; tg® x, + 


Beta ee) 
— |B, Ar Sin (tg x’) — const] cos x’ 


/ 4 
1 = (2B, tg x’ + = Batg? x’ 4 


6 8 
+ 5 B, tg® x’ + > Bs tg? x’ + =) 


(2,9) apes 5 
+ 7 B, tg? x! alle E Be tg oll aE 


7 
+ 6 B, tg® x! -+ «| 
+ [B, %Xr Sin (tg x’) — const] sin x’ 
worin noch — wir bekannt — 


Mr Sin (2 7") ==Aeein at = 


Y 
u=0 ul (2a 4) yrs 


Abbildungen 3, 4 und 5 geben eine Darstellung der 
Strahlen und Wellenflache eines komabehafteten 
Strahlenbiindels fiir 3 verschiedene Falle. 

Untersuchen wir noch die Linge der Lichtwege eines 
beliebigen der Strahlen der Meridianebene und des 
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Schwerstrahls (Hauptstrahls) vom Objekpunkt bis 
zu ihrem gemeinsamen Schnittpunkt im Bildraum. 
Zu dem Zwecke greifen wir wieder einen beliebigen, 
dann aber bestimmten der Strahlen heraus. Sein 
Neigungswinkel gegen den Hauptstrahl sei x. Auf 
ihm tragen wir vom Schnittpunkt mit dem Haupt- 
strahl die zugehorige Lange — A di ab und betrachten 
die durch den Endpunkt der abgetragenen Strecke 
— Al; gehende Wellenfliche bzw. deren Schnittkurve 
mit der Meridianebene. Fiir diese gilt 
Bilao 
sin x, 


(2,10) 


so dass fiir diese Wellenflache die Konstante nach 
(2,9) den Wert 


CE) Cohn = Al, -- 
(FB wae 4 6 ; 
aa 2B, 4 | B, tg? x, t : B, tgt x, i 
cosx, \ 3 5 / 
| ) y 5 
[eB tos, 4 By (eee 
COS %, 2 


+ B, Ar Sin (tg x1) 


Die Gleichungen dieser betreffenden Wellenflache 
lauten mit dieser Konstanten (2,11) nach (2,8) bzw. 
(2,9) demnach 


ph Bh Sees 
\ COS x, 
py es 4 cos x! I 
4 Gage ie tat) + 
‘ COS x, 
j 6 cos x! | 
-. B, tot x! =F oa Bs as tg* * Si 
5 COS %, / 
a 3 cos x! 1 
+ | By tee ee ae tg x, \+- 
\ 2 COS X%, 
5 COS tte, aorl 
+ ( B, tg® x’ + — B,; = te? oc) apes 
4 COS x, 
+[—B, ( Ar Sin (tg x’) —%Ar_Sin (tg x) ss i cos x’ ; 


Lie E B, [tg series = + 
\ COS x, | 
/ x! 
> — Bi | te?) ; tg ay 
i COS %, 
ii if 
+ — Bg / tg® x’ — tata ae 
COS %, 
3 sin x’ ' 
+ | — B,{ tg? x’ — , tg %, 
} 2 COs %, 
(eae 
A COS %, 


+ [B, (x Sin (tg x')—Ar Sin (tg x.) ) + Al] sin x’, 


\ 
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Ass. 3. — Strahlenaberration, Strahlenverlauf und Wellenflachen [im Schnitt 
mit der Meridianebene des bildseitigen Hauptstrahls ( = Schwerstrahls) ] 
eines mit Koma behafteten Strahlenbiindels fiir den Fall, dass 


1 
INE es ae Signe & OF = = baw. 1 000. 
) 


Fir x' = 0 wird 4, =(4-:)»9 = 9 — wie erfoderlich— und Setzen wir hier wieder den Wert fiir AL aus \(2,11) ein, 
so ergibt sich fiir x’ = 0, also fiir den Schnittpunkt der 


2 Be Bekele, I \ ; ! 2.' .  betreffenden, durch x, und AL, bestimmten Wellen- 
fr = (1)o = — 2 2 | = Aelia SR Hew OK . ‘ 
x COS x, / cos x, flache mit dem Hauptstrahl des Biindels 
6 pad (Eo = — 2B, + const,. 
-B,+ —B,tg?x, + >1/0 See 1 
8 Fiir einen dem Schwerstrahl (Hauptstrahl) unmit- 
=| 7 By tg* x, + | telbar benachbarten Strahl gilt x, > 0 und Al, > (0, so 


dass fiir eimen solchen Strahl die Konstante (consty) 


l lits ° der nach der gleichen Vorschrift wie oben konstru- 
. er | ° B. + 2 Be te2 x’ + der nach der gleichen Vorschrift wie oben konstru 
cos x B% Tad 7 Fed 8 ierten Wellenflache nach (2,11) den Wert 
a! li sonst, = + 2B, 
= & Be tg? at = =| : BG : 2 
' annimmt, also (&9)) = 0 wird, der Ursprungspunkt des 


+ By %r Sin (tg x,) — Ad,. (€, e) — Systems also im Schnittpunkt des Haupt_ 


1. PICHT [ Orr. Acta 


| 0.6 


O02 


+l’ ~~ 1 & ( 
Be SPS -—| = Es +& 
ae MS ews a 
PON 
| -O.2 ff tS 
--04 / 
L-06 
Ass. 4. — Strahlenaberration, Strahlenverlauf und Wellenflachen (im Schnitt mit der Meridianebene 
des bildseitigen Hauptstrahls) eines mit Koma behafteten Strahlenbtindels fiir den Fall, dass 
@) 
ANE on ieee! See Ee! ly z bzw. 1000; ag = — — bzw. — 4 000. 
) vd 
strahls des komabehafteten Strahlenbindels mit den . : ' ewes 4 ! 
ihm unendlich benachbarten Strahlen liegt. eee Wa ysk | pea Si 3 Bites 
Ware const, = const), so wiirde es sich bei beiden Fa 
Wellenflachen um die gleiche Wellenfliche handeln, 6 ; 
der zum Hauptstrahl geneigte Strahl (mit Al’ = Al,) a = Be tgt x +... 
sowie der Hauptstrahl selbst hitten dann also in y 
ihrem gemeinschaftlichen Schnittpunkt gleichen Licht- ms | ee te 5c =) B, te? x) 4+ 
wegabstand vom gemeinschaftlichen Objektpunkt. cos x, | 2 S SRT a 
Tatsachlich aber ist const, ~ consty. f 
: : e , . + B, % Sin (te x,). 
Der Wertunterschied beider Konstanten gibt wieder ia (8 4) 
unmittelbar den Lichtwegunterschied beider Strahlen . : 
von dem beiden gemeinsamen Objektpunkt bis zu AWijo = | - : st) QBs Ali) 
ihrem gemeinschaftlichen Schnittpunkt im Bildraum. \ cos x, 
Wir bezeichnen diesen Wegunterschied wieder durch ' 
/ ye Bae : B Batt i oll 
AWajo = W (tg x.) — Wo = (Ero — (Lolo = (Eile ggg ag gt eet Oss 0 
1 . 


Dann ist nach (2,11) + B, tg* x “five )ape oy tek Sm Go %) . 
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¢ bgk’ 


Heuptstrah #é 


Ass. 5, — Strahlenaberration, Strahlenverlauf und Wellenflachen (im Schnitt mit der Meridianebene des bildseitigen 
Hauptstrahls) eines mit Koma behafteten Strablenbiindels fiir den Fall, dass 


1 2 
Al’ =a, tg x’ + a, tg*x' + a, tg*x’ ; a, = Ge bzw. 1000 ; a, = —— > baw. — 4 000; a,=1 bzw. 6 000. 


Schlussbemerkung. — Bei den vorstehenden Unter- 
suchungen, die fiir die Beurteilung der Korrektion LITERATUR 
eines optischen Systems von Bedeutung sind [2], [1] J. Prcnt, ZS. f. Instrkde 51, 1931, 77-81. 
handelt es sich allein um die Bestimmung der Licht- : Ra : 

2 aoe ; 2 [2] Siehe das inzwischen (Dezember 1955) im Verlag der 
weg unterschiede. Dieindem Gebiet der Brennflachen Wissenschaften, Berlin erschienene Buch des Verfassers : 
bzw. beim Durchgang durch die beiden Kaustikschalen « Grundlagen der geometrisch-optischen Abbildung » 

— von denen die eine evtl. entartet ist — auftreten- Manuscrit recu le 12 novembre 1955. 
den Phasenspriinge sind hier nicht beriicksichtigt. 
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A procedure for the differential correction 


of optical systems allowing large parameter changes. 


W.S. S.. BLASCHKE. i 
Cooke, Troughton & Simms, Ltd., York. * 


Summ+ry. — A method of differential correction for opticul systems is developed which serves in the coarser analysis of the aberrational 
characteristics of the system as well asin the final detailed balance of residuals. © 

Variations of optical paths through a general oplical system due to relatively large Gee Poem changes are, considered, both 

al the changed refracting surfaces and surfaces subsequent to them. The extra-axial meridional case is tceated as well as the “axial 


case. 


t 


The technique is applied to two high numericgl aperture microscope objectives suffering from the marginal ana zonal unbalance 
normally associated with such objectives. A substantial change in the mranetnont shapes is achieved, resulting in y much improved 


correction. 


2 correction dijférentielle des systémes optiques utilisable pour une analyse grossiére des aber- 


rations dun ae ausst dien que pour la‘correction finale des résidus daberrations. 

Les variations de chemin optique a travers un sysléme ayant subi des modifications relativement importantes, dés paramétres 
de construction sont éludiées a la fois sur les surfdees réfringentes modifiées et sur les surfaces suivantes. On éuldie dans le plan 
méridien, le cas dun point objet éloigné de Vaxe aussi bien que celui de Vobjet situé sur Paxe. 

La méthode est appliquée a deux objectifs de microscope de grande ouverture. numérique, affectés de défauts d’équilibre 


des aberrations marginales et zonales propres @ ces sorles d’objectifs._ 


est calculé ; il améliore sensiblement la correction. 


Un changement: nolable des formes an front donde 


# 


clion oplischer Systeme entwickelt, die sich bei der Probena- 


aie des areca bie gut bewdhrt wie bei der Feinkorrektion der Réstfehler. Die Anderungen der optischen Wegldnge durch 
ein allgemeines optisches System bei verhdltnismdssig starken Anderungenim Aufbau, werden untersucht und zwar sowohl an der 
gednderten brechenden Fldche als auch an den anschliessenden Fldchen, Der Meridionalstrahl wird dabei ebenso gut erfasst wie 
der achsiale Strahl. Dieses Verfahren wird auf zwei Mikroskopobjektive hoher numerischer Apertur angewendet, bei denen Rand 


und Zone, wie bei solchen Objektiven iiblich, nicht ganz ausgeglichen sind. Man erreicht eine wesentliche Anderung in der Form 
der Wellenfront und damit einen sehr viel besseren Korrektionszustand. \ 


I. Introduction. — The fine correction of optical 
systems based on calculations of increments in aberra- 
tions with small constructional changes is today a well 
established procedure of optical design. The technique 
has been presented in various forms (CrRuIKsHANK & 
Mc AULEY, 
DERT, 1954: 
the latter three authors for automatic computation. 
Throughout all these papers the difficulty presented 
by the non-linear dependance of the aberration on 
the parameters and the limitation of the USOT ED to 
very small changes is accepted. 

The problem has been stated recently in a different 
way by Hopxrns (unpublished) as follows. If a small 
change is made in the construction of the system, say 
a change in the curvature of a surface, the wavefront 
incident upon this surface will travel in the modified 
system, at a particular zone, a greater distance in one 
medium, and a correspondingly smaller distance in 
the adjacent medium, than before. This introduces a 
change in the departure of the wavefront from a refe- 
rence sphere and hence a change in aberration. If 
however the parameter change is a first order small 
quantity, the path taken by the ray in the modified 
system will be neighbouring to that taken by the 
original ray. Ray paths will be equal by FEermat’s 
Theorem and the change in wavefront departure from 


* Now at Taylor, Taylor & Hobson, Ltd, Leicester. 


1945-46 ; Freitac, 1952; Rosen & EL- . 
Lucy, 1955), and has been adapted, by 


reference sphere may be transferred to the exit pupil, 
The argument is completed by noting that the paraxial 
focus has moved and that a change in reference sphere 
to the new focus is needed to obtain the required value 


of the aberration change. This approach leads to much 
simpler equations than.those in the above papers, 


while making exactly the same approximations. 
Moreover it lays bare more precisely the nature of 
approximations made and shows how the technique 
may be modified to give the aberrational increments 
for very much larger constructional changes. It is the 
purpose of the present paper to develop this extension 
and thus provide a method which will assist in the 
coarser analysis of an optical pr as well as in 
its final correction. 


IJ. Treatment of Approximations. Let = be an 
optical system with ray OTI (wavefront 6 I) incident 
upon it. At some point inside the system we have a 
small constructional change, e. g. an increase in the 
curvature of the refracting surface T; which becomes 
I's. Let N; and N; be the refractive indices of media 
on the object and image side of I; respectively. 
The ray path OIIP;1’O’ becomes OIP;P*II’*O'*. II’ 
and Il’* are the points of incidence of the original and 
varied rays at the exit pupil reference sphere whose 
centre is at o’, the original paraxial focus. P*Q; 
the perpendicular dropped from Px to the ray P, ...... II’. 


vol. 3, n° 1, March 1956] 


If by the theorems of Matus and Frrmar the 
paths, P&+...5..04. TRIE RC as beg ene II’ are equal, the 
change in wavefront aberration taken as the depar- 
ture of the wavefront from the reference sphere cen- 
tred at o’ is simply — [NV,(PP*) — N‘{(PQ)]. The 
reversal of sign is in accordance with the convention 
that a longer path along the ray, which means a 
wavefront retarded with respect to a fixed reference 
sphere, gives negative-sign aberration. Reference to 
o’*, the new Gaussian focus, then requires adding the 
extra path N’o travelled along the ray between the 
reference spheres centred at o’* and o’ respectively 
(N’ is the refractive index of the image space medium). 

There are two basic approximations here which 
break down to a greater or lesser extent according to 
the position in the optical system of the constructional 
change, its magnitude and nature, and according to 
the position of the exit pupil. The first approximation 
is the assumption of neighbouring paths already men- 
tioned. The second assumption is that the original and 
varied rays at the exit pupil are so near that the change 
of focus term (path between reference spheres centred 
at o’ and o’*) may be measured indifferently along 
either ray. 

The first approximation becomes much less severe 
when a point is found in image space on the original 
ray which has the same optical distance from Ba and 
Qi to an accuracy which includes the second power of 
the angle between the rays. Such a point is found by 
tracing the paraxial tangential fan about the original 
ray as axis, which marks the progress of the wavefront 
element P* Q; of aperture w,/;, now considered no 
longer a first order small quantity. The trace of such 
a paraxial fan commonly used in astigmatism calcula- 
tion involves very much less calculation than the 
retrace of a trigonometric ray through the changed 
system. We are at the same time led to the elimination 
of the second assumption. The new focus O’* is to the 
same approximation the point where the tangential 
fan ray intersects the axis. If we drop all reference to 
Gaussian foci as such, we can easily refer the aberra- 
tion change, given with respect to a reference sphere 
centred at O’*, to a sphere centred at O’, the original 
ray focus. The latter focus has been called (HopKrns, 
1952), the invariant focus, and has the advantage of 
giving a unique value for the wave front aberration, 
ij. e. as the wavefront progresses along the ray. its 
departure from a reference sphere centred at O' remains 
invariant. All reference to the exit pupil of the system 
is thus avoided and the second approximation is 
circumvented. 

The technique of correction is applied as follows. A 
sufficient number of trigonometric rays, axial and 
extra-axial, are traced, say 15 : marginal, sub-marginal 
and zonal on axis and for 2 positive and 2 negative 
object heights, (the latter are traced as « T-rays », 
(see below) relative to trigonometric pupil rays). The 
correction of the system is satisfactory when the Opti- 
cal Path Differences, axial — ray path, (henceforth O. 
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P. D.) at these 15 foci are simultaneously zero. Diffe- 
rential changes are calculated on each of these rays 
and so combined as to minimize the 15 quantities. If 
large changes are involved which modify the basic 
structure of the system or give rise to completely 
different ray paths, it will be necessary at a certain 
stage to retrace the system and repeat the process. 
Simultaneously zero O. P. D.’s at all invariant foci 
would result is spherical wavefronts centred at paraxial 
and paraxial-tangential foci. Thus for a satisfactory 
result, attention must be given in the earlier stages of 
the design to the requirement that the system should 
have as far as possible a flat tangential field (as given 
by conventional tangential fans). 

We shall develop expressions for the aberration 
increments working for simplicity first for the axial 
case, showing later how a small number of auxiliary 
calculations make them adaptable to the extra-axial 
case. First we must consider the change in O. P. D. 
at a surface when the ray path through it is varied. 


III. Effect on aberration at subsequent refracting 
surface of a constructional change. — I is a refracting 
surface (pole A) between media of indices V, NV’. PO, 
after refraction (henceforth a. r.) PO’, is a ray 
through the surface, P being its point of incidence. 
Reference spheres are drawn through P with O and 
O' as centres intersecting the axis at R and R’ respec- 
tively. AR = «, AR’ = e¢’, with the sign convention 
that ¢(c’) is negative when at the axial intersections 
the refracting surface is in advance of the reference 
sphere and vice versa. Then: 

O.P.D.o....0' = N(AO) — N'(AO’) — N(PO) + 

+ N'(PQ’) = — N(AR) + NAR’) = A(WVe), 
where A denotes the difference (expression a. r.) — 
(expression b. r.). 

We now take the ray PO, a. r, oBe to be that 
traced through the unmodified optical system. A cons- 
tructional change has been taken at surface Ij (fig.1) 
before the surface I. of figure 2 P*Q; and ut; have 
been obtained there and a tangential fan has been 
traced with these quantities as starting data. The ray 
from this fan is (in fig. 2) P*O* a. r. P*O’*. It deter- 
mines with the trigonometric ray foci 04, of ; apertures 
u,, Ww, and incidence heights y;, y{. We call 

Po St, Pol St 3 POSS , POEs"; 
HOLS 5 LORS U 


and the perpendiculars from o, and 0; to the axis 
respectively ¢,, et. \ 

We have equality of optical paths between o, and 
o; along the two rays and draw reference spheres 
through 0; and 0} with centres O, O* and O’, O’* respec- 
tively. The spheres determine axial intercepts BB* 
angebro*: 

Evidently the change in O. P. D., O%*...... O'* (via 
exis)——-viat) trom.O, P: D. O....;. O' (via axis) — 
(via P) is 
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(1 a) A(Ne)* —A(NVe)=—A (ys a aS at 


nl 4 sin U 7 
} heya aie oeaey tere 


=— A(M(S — t) ® (sin VU, u,)) 
(1 a) remains finite in the two cases : 


(1) A(Ne)*—A(Ne) = A [V(BB*)] 


\ 


‘sin U+sin u, 
=4(v)s—1—(s—a| 1) 


sin (U + ut) 


This becomes after some transformations which 
neglect powers of u; greater than the second (see also 
Appendix I) Image at infinity : S’ > o, sin U’ = 0, as 


Fia, 1. 
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/ Ve but intersecting in P* instead of oj. The corresponding 
1b N'ce')* — LAIN 5 area eh T t 
) oe es) g oe Mt axial intersections are R* and R. We have 
e* — ¢ = PP* oe AR*, (ce =A =O: 


? 
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and tangential focus at infinity : t’ > 0, ul = 0as 


(1c) (N’e')* — (N'e!) = N' a ; Therefore 
ace (3) A(Me)* —A(Ve) = N'(PQ) —N(PP*) — W'(RR*) 
The expression ({ a) however does not have to be Q is again the foot of the perpendicular from P* to 
computed at every surface in the system since PO’, which may be confused with the portion of the 
pad i EP ee TG reference sphere centre O’ through P*. Q* is the inter- 
v's) ee UNS) Suey section of reference sphere centre O’* with PO’. 
It follows from what has been said about invariant 


where j and j + 1 denote consecutive surfaces. In 
fact after the calculation has been completed at the foci that N’(RR*) is equal to the term N’(BB*) of (4) 
surface whose constructional parameters are varied, § III, when it is remembered that 
say surface i we have for the complete result, Q* o! 
P* of = + O(ut). 


(2) 3}OPD) = S(OPD); + [N(S — t) sin UV, us)li+4 Qot 
—[N' (S’ —1t’') ® (sin U'," u! 
re ) Peas Ya ahs Thus the last term of (3) cancels with the second 


where k denotes the last surface. The difference of the term of (2) § III. We proceed to derive an expression 
last two terms is the quantity assumed vanishingly for the first two terms of (3) and for the y; and u; 
small in the approximation of Frermar’s Theorem which mark the tangential fan to be traced through 


(§§ I, II). It is the complete expression for the O. P. D. the remainder of the system. 
change subsequent to the modified surface. In what follows rigid adherence to the same degree 
of approximation has been abandoned. Thus terms 


It remains to find expressions for the parameter 
derivatives of the O.P.D. at the surface which is depending on the square of the parameter change are 
modified. The cases to be considered are: (1) achange neglected where these can be considered random (}). 
in surface curvature and “‘ bend ” of a component ; There is however one case where squared terms are 
(2) a change in position of a surface or component; of particular importance i. e. in the ‘‘ bend” of a com- 
(3) a change in refractive index of a medium. ponent near minimum spherical aberration, the vertex 

of the parabola (Rumsry, 1951). The non-random 
square terms are then separated out and independently 


IV. Change in curvature of refracting surface and 


« bend » of component. — In figure 3 we consider the 
refracting surface 1 of figure 1 (§ I). Reference spheres 1) Apart from simplification of formulae there is an advan- 
are drawn as in figure 2 with O’ and O’* as centres tage in the results being directly scalable. 


computed. 
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Using the formulae of the ray trace scheme given 
below (Appendix II) we have: 

— 2hsinl 
GF=> G 5 


which gives upon differentation with respect to c, after 
rearrangement, 


G = (4 + cos U) (14-+c0s 7) + sin U sind 


2h?8e (1+ cos U)(1+ cos l) 
GeosI ° G 


(4) PP*¥=e*—e= 


but PO: = PP* cos(1 — J’), and with the aid of 
A(N sin J) = 0 

(5) N’(PQ)— N(PP*) = (e*—«) cos J A(N cos £), 
while 
se A(N cos I). 


Let /* (b. r.), //*(a. r.) be the angles of incidence 
of OP* a. r. P*O’* with '*. As U is unchanged, w; is 
given by the change in the deviation of the incoming 
ray. Therefore u, = — A(J)* + AW) = — 8} A(/)}. 

But in the form of ray trace given in Appendix II 

sin l = he — sin U 


(6) y! = PP* sin(J — I’) =(c* 3S 


where f is the perpendicular from A to OP. 
Therefore 


h&e 
7 oi =O) 
“) COS7/ aan 
while because of A(W’sin I) = 0, 31° = 08 
N' cos 
Therefore 
N cos I 
ba Oy SEE eee 
(8) mt ( N’ cos I’ ) 


(8) may be interpreted as follows: The change I, I’* 


means at the point of incidence of the ray OP a change 
in the tilt of refracting surface relative to incoming 
wavefront. If we tilt instead the wavefront by 2=—<8J; 
the new ray is treated by T as was OP by I™*. 
Thus a fan with apertureo = Q, y, = 0 is traced 
through the surface and the resulting w’ is referred to 
the ray PO’ by subtracting from it the tilt angle ®. 
If the subsequent surfaces are unchanged we use mt 
and y'{, of (6) to trace a fan through to the end of the 
system (?). ; 

This interpretation of (8) leads directly to a method 
of calculating the ‘‘ bend ” of the component. Such a 
“bend ”’, we take not as equal increments 8c imposed 
on all surfaces but as a constant tilt 2 of the incoming 
ray. The effective changes of curvature are then given 
for each surface by 

2 cos I 
(7 a) h 


Such a procedure is known from first order calcula- 
tions based on a paraxial ray. However, while in the 
latter case a constant increment ydc (y = paraxial 
incidence height) gives exactly unchanged deviation 
of the paraxial ray this is only approximately so in our 
case. Further, while in the paraxial case the increment 
ysc may be scaled for other rays in the finite case Q 
can be constant for only one ray and has to be obtained 
anew from (7) for the other rays. 

In figure 4a we consider two consecutive surfaces 
of a component. They are “ bent” by Q,, and Q, cor- 
responding to Sc, and Se, [by (7)]. The 8c, and 8, 


sc = — 


(?) In working out TE 8 in Ut — Du; we are justified in 
the use of D, the distance along PO’ to the next surface in 
Ryne of D — PQ since only a square term is thereby neglec- 

ed. 
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are in turn derived from a constant [by (7a)] at 
some other ray. The path 1, 2 through the unbent 
component becomes 1*, 2* through the bent component. 
We require the difference i in path from P,, via 1, 2 and 
via 1*, 2* to a point 0, which is found from a ‘‘ tan- 
geritial fan ” traced in the Loans manner. 

a) A tangential fan with o, = Q, ; yu = 0 is 
traced aes A 

b) i = wo, — Q,. 

C) Yr, is obtained from (6). 

d) ee Yt, — Di? Uy! where D, , = P,P,. 


e) Oo = = Uy a Qo. ! 
f) ‘A tangential fan with », y;, 18 traced through 2. 


I 


g) Uy, = @, — Q». 


¢ . : cos I} 
h) Vig is found by adding to the value y, ———— 

2 cos I, 
derived from the tangential fan in (f) the result of 
applying equation (6) at surface 2. 

i) If the component subjected to ‘‘ bending ” ends 
at surface 2, we have the required point oj (in this 
case 0,’) at t’ from surface 2, where y,’ = t/ u,! is 
the value obtained in ‘h). t’, Uy are then directly 
used in equation (2) of § III. If ‘‘ unbent ” components 
or surfaces follow, a conventional fan is traced on to 
the end of the system and (2) of § III is applied in 
image space. 

j) If more surfaces remain in the component subjec- 
ted to ‘ bending ” or if the following component is 
also “‘ bent ’ the procedure (e), (f), (g) and (h), for 
such surfaces is repeated. 

Refering back to figure 4a we divide path difference 
(N. B. sign convention in accordance with § II). 
#2 — 1%, 2* into components 1,°2 — 1*, 2° and 
1*,2 — 1*, 2*, that is, we consider the change in cur- 
vature of each surface in turn. The procedure outlined 
above in steps (a) to (j) is to relate to this construction. 


We have 


ee oe pt Le Pl 
Ves) VPP *) 

and 
4¥3£1,2* —[p, ——— = Q¥]—[p,-———. PH] = 
= N3(P2q5) == (P ps) + Of de?) 

but 

y= — Q*] = (03 ==> al: 
Hence 

V2— AF 2* = [Py q.] — [Pi = Pre 
= »; (<* — e) cos J A (N cos J). 


The point o,' is equidistant from Py and q, and 
P¥q, is the yy, obtained in step (h), since 


P¥q, = poQ, + Pray + O(8c?). 
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Also ~ P¥ o,! q, is the required angle uw‘, except for 
a term depending on 8c?, being the addition of the 
angles made by path 1*,2 with 1,2 and 1*,2* with 
1*,2. Here however step (f) does not correspond 
exactly to the construction of figure 4a viz. the fan 
marking path 1*,2 relative to path 1, 2 through sur- 
face 2 should be traced with starting data Ytyy Ut, not 
yt,, @ asin (f). The additional path difference involved 
can be seen from figure 4b to be a quantity depen- 
dant on 8c?. OPO’, figure 4b, is a trigonometric ray 
through a surface. P*o,, a. r. P*o; is a fan with 
apertures w,, a.r. uj. of P*, a. r. 0,*P* is a fan with 
aperture @,a.r.@'. o = u, — Q, wo = u, — Q. We 
require the change in the path of o* to o/* from that 
of 0, to of via P — via P* in each case. The conside- 
ration is analogous to that of § III and leads to the 
form (1b) § III with y, taking the place of p,and Q 
the place of uy. 


1 4 
Bae QA(Ny,) = =e A(N cos 1), 
where z = i = i ,and has to be subtracted 
cos [ cos I! 


from the path difference 1,2 — 1*,2*. Importance is 
attached to it as a non-random square term. One of 
the characteristics of a ‘‘ bend” is the progressive rise 
or drop in incidence heights of a ray through the 
component according to the sign of 8c. Thus QzA 
(Neos J) will never average out. It becomes in fact 
the controlling term at the minimum of the spherical 
aberration parabola where 0(O. P. D.)/22 = 0. Sum- 
ming up, we have for a “‘ bend ” of a component, when 
the calculations are carried through as above, 


(9) A(Ne*) — A == 
—e)cos!A(NcosJ) + Es ZQA(N cos I) — 


V. Displaced component or surface. — The pro- 
blem of an axial displacement of a component or 
surface is very simply treated by considering a con- 
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jugate shift in the opposite sense. Two such conjugate 
changes will reproduce the axial displacement of any 
part of the system. What we are deriving is an exten- 
ded form of Herschel’s condition in which displace- 
ments are considered significant to the second power. 
OR P’O’ is the original ray through the system 
Y. Let the system be axially displaced toward O by 
an amount 8z. The corresponding conjugate change is 
— 8z. A fan traced with u, = 0, y, = — 8zsin U repre- 
sents this accurately. If uj, and t, = P04, are 
determined by this fan in image space we have 


(10) A(Ne)* —A(Ne) = 
= — A[N(S —t?) Osin U, u,)] 

U 
=—Ny; me 7 —[N'(S! 1") (sin U’, ui), 


Where on account of u, = 0 the form (4 c) of § III has 
been used in object space. 


VI. Change in refractive index of medium. — In 
figure 6 we have surfaces 1 and 2 of a component 
surrounding a glass of medium NV, = JN, whose 
refractive index is changed by 8N, = 8N, = SN. 
Using a symbolism similar to that of figure 4a we 
have a path 1, 2 through the unchanged system ; 
path 1*, 2 through surface 1 between media N, 
and Ni + 8; path 1, 2* through surface 1 between 
N, and Ni, andl total surface 2 between V, + 8N 
and NV}; path 1*, 2* through the lens of changed refrac- 
tive index. 

pP*O;* (fig. 6) is the final path 1*, 2* and intersects 


PO! (pat 1, 2) in 0}, determining Ploy = t, and 
angle u,’. We find this data as follows : 

a) The angle u,' (angle between 1*, 2 and 1, 2 after 
surface 1) is found from 
We d(A(Z)) = 1, — 1* 


where Ni sin I, = (N, + 8N) sin I* 


and therefore 


SN tan J, 
(11) uy, = —— 
N, 
b) 
Du yi; 
(12 yeenn 25 2P iS ee 
cos/, cos J, 


A tangential fan is traced with es and z, through 
surface 2 between media NV, and N, Gath 1*, 2) giving 
angle Waa 

Cc) Uta the angle between paths 1, 2* and 1, 2 after 


surface 2 is found from 
Wop = I,— I;*, where (N, + 8N) sinJ, = 
while A(V sin /,) = 0 and therefore 


Uf bd I 
JV, sin pe 


8N sin J 
(13) Ui ah eee 
NV, cos I, 
a) —* where 2p ut Ui: A square term 


ty 
arising from differentiation of the tangential fan, b) viz 


Fic. 6, 
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N,8N 
su, = ——__ 
za N, cos I, 
significant. The difference in paths 1, 2 and 1*, 2* 


(respecting sign convention, viz.sV > 0 means a longer 
ray path) then becomes 


A(N,¢,)* oe A(N, €,) + A(N,e,)* a A(N, f,) = 
= SN(c, — ¢,) +N, 8c, -- N, Se, — N, 8e, 


(wi... + 2,¢,) may be added if 


where the Se values are the axial intercepts marked 
in figure 6. But by the property of invariant foci 
(O,, O,*) 32, = 3c, while N, 3<, again cancels with 
the second term of (2) § III and «/ — «, is equal 
to d-D, 1. e. axial separation of surfaces minus separa- 
tion along the ray. The complete result reduces to 


(14) A(Ne)* —A(Ne) = 
=e = 1) ONAN (1S! —t) @ (sin U7) 
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on the pupil ray, the O. P. D. is again given by A(Ne). 

In figure 7 we show a typical pupil ray, T-ray 
through a surface and the variation of both rays. The 
constructions and notations are obvious extensions 
of the axial case. 

If we now repeated the considerations of § 3 we 
should arrive at an expression analogous to V(S — t) 
(sin U,w,) in the form 


(15) N [(S — 4) D(sin Vem) =- 


Ue 


i =scosma 
(oS (eae Py) coma | 


+ (L—t?) O(sinV, u,)— 


The proof of (15) is laborious and is not given here. 
Instead both the formulae are treated again in Appen- 
dix I by a more elegant method. So far as the considera- 


Fie. 7. 


VII. The extra-axial meridional case. — It is 
conventional practice to treat the extra-axial (meri- 
dional) aberrations by tracing rays through the point 
of the object field under consideration which are in 
effect axial rays imaging fictitious object points. For- 
mulae for the extra-axial O. P. D. have been given by 
Conrapy (see Hopkins, 1950) which introduce data 
from such rays together with data from the pupil 
ray (which images the entrance pupil). As in off-axis 
imagery one is only concerned with the relative paths 
along the two rays, it is reasonable that the extra- 
axial ray should be traced in the first place relative to 
the pupil ray (T-ray) and not the axis. (Formulae 
for such a ray trace scheme are given in Appendix If). 
If such a trace includes quantities ¢, <’ defined as before, 
§ 2, but with the axial intercepts becoming intercepts 
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tions of §§ IV, V, VI, are concerned the simplest 
procedure is to extract from the traced T-rays, the 
quantities used in these sections, which are then consi- 
dered relative to the axis as before. Thus after the 


obvious derivation of sin U, cos U from U = U + V 
we use h = (sin U + sin J)/c, all other quantities 
appearing explicitly in the T-trace. The pupil ray 
itself is an axial ray and its variation requires no addi- 
tional considerations. 


§ VIII. Application and discussion of numerical 
results. — The system of differential correction has 
been applied to two designs of microscope objectives : 
A 2m.m. Oil Immersion Achromat of 1.30 N. A., anda 
4mm dry Achromat of. 85 N. A, 

In interpreting Tables I and II we should note the 


iss 


[Orr. Acta 


its) Ww. S. S. BLASCHKE 
SL 72-362 ins §Z- .609ins SL/= .259ins SL%.05& 
Axis MARGINAL !1.30NA Axis 
T= -0035 ins | ae 2 
3.09 -046 SUB-MARGINAL 326 225 


H>9 


| : 
5 498 ZONAL ! 


we =i 


Original 
system 


2 R) 4 
OPD.9/(WWavelength) 
Modified 


system 


Y=-.0095 


Vignetted at\/..710 
objective front 


IN 1.04! -6.68% 
ee 
SL5--5.12 
Fie. 8. 


following : The O. P. D.’s at invariant foci are always 
opposite in sign to those at the paraxial focus. The rea- 
son for this is that the change of focus term predomi- 
nates over the path difference along the ray. Thus 
an undercorrection (positive at the paraxial focus) 
is indicated by a negative O. P. D. at the invariant 
focus and vice-versa. Comparing at a given aperture, 
the O. P. D. on axis with those derived from the T- 
rays, we see that the values for negative and positive 
fields (denoted by the sign of H, the LAGRANGE inva- 
riant) lie on either side of the axial value. This indi- 
cates an odd term in the extra-axial O. P. D. which 
expresses the coma. On the other hand the non-paraxial 
astigmatism and excess of spherical aberration on the 
principal ray over that on the axis are expressed by 
the mean of the extra-axial O. P. D.’s. (at invariant 
foci) being greater (overcorrection) or smaller (under- 
correction) than the corresponding axial value. 

In Table I results are quoted for seven construc- 
tional changes of the first objective. They are taken 
in two stages and so chosen as to leave the convergence 
of the marginal ray unchanged (i. e. Buy ==.) ocA 


comparison is made by retracing the system after the 
first three changes. The remaining four changes are 
calculated on the intermediate state and compared 
with a retrace of the objective in its final form. 


In its original state (see also fig. 8 a) the objective 
shows two defects : 1) the axial wavefront, while very 
close to the paraxial reference sphere (in fig. 8 a the 
y-axis), deviates rapidly at the margin viz: has a 
rapid change in gradient ; 2) the coma is overcorrected 
or outward (negative in our convention i. e. wavefront 
for H > 0 retarded with respect to reference sphere 
while in advance for H <0). In Table I, 1) is seen from 
the large negative O. P. D. at the marginal invariant 
focus ; 2) from the fact that the O. P. D. at the sub- 
marginal T-ray focus is more positive for H > 0 than 
for H < 0. 


The result of the modification has been to eliminate 
the marginal under-correction at the sacrifice of a 
small amount of zonal undercorrection. The latter is 
largely taken care of by the choice of best focus. Thus 
in figure 8 a the wavefronts (referred to paraxial foci) 
before and after modification are plotted against the 
square of the aperture. The paraxial reference sphere 
(also the reference sphere for all foci lying in the Gaus- 
sian image plane) is taken as coincident with the y-axis. 
The reference sphere at any other focus is a straight 
line inclined to it. Hence the asphericity of the wave- 
front is judged by its approximation to a straight 
line. 


The coma has been largely removed. The O. P. D.’s 
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TABLE I 
. . me 
Objective : 2 mm. 1.30 N. A. Axis Fields 
Jef Sony vai 0 
Achromatic Lens. Mgl. S. Mgl. Zonal =e e 
Mgl. S. Mgl. Mgl. S. Mgl. 
Original System O. P. D. — 00 0094 00 0004 00 0003 00 0191 00 0144 — 000051 00 0003 
DO Nee an 2 siete ds «eve oo tha Fc 0577 0216 0040 0949 0426 0337 0101 
Xs ==> pW SSE ae Genon moe — 0420 — 0202 — 0047 — 0692 — 0320 — 0226 — 0114 
C14 BOS Or ateseacus te avods homer 0005 0003 0001 0003 0002 0024 0021 
3(OPD) by differential Correction ..... 0162 0017 — (0006 0260 0108 0135 0008 
Intermediate System a 
OPD* (by differential Correction) .... 0068 0021 — 0003 0451 0252 0084 0011 
Pr) Sa(TECLACEC)) ae miewsait «ei siere oes as 0007 0009 — 0001 0191 0128 0168 0064 
Pet Og COTTA Ta ah ein oon does 0314 0119 0023 0476 0184 0233 0078 
dc, Sas ASS euciche tenia aiion stayo vee ia — 0311 — 0132 — 0030 — 0527 — 0212 — 0177 — 0076 
025-9 == Qe n a tapenah toms, «, cee oy sae, 2 — 0003 — 0003 — 0003 — 0079 — 0078 0 — 0009 
8Cy4 OT 9OLS Rererece rs seis + Fetes 0012 0009 0003 — 0022 — 0018 0032 0026 
S(OPD) by differential Correction.... 0012 — 0007 — 0007 — 0152 — 0124 0088 0019 
Modified System 
OPD* (by differential Correction) .... 0019 0002 — 0008 0039 0004 0256 0083 
ONPD K+ (GRABLE GON in S ae slo cabin A oI ea cin 0002 — 0003 — 0012 0073 0057 0233 0089 
TABLE II 
Objective: 4 mm. 85 N., A. Axis Fields 
Ha) Jobe MY 
Achromatic Lens. Mgl S. Mgl. Zonal 
Mgl. S. Mgl. Mgl. S. Mgl. 
Original System O. P. D. 00 0128 00 0066 — 00 0048 00 1026 00 0346 — 000536 | — 000177 
O21 Se (VOGT Die st eyes,s/okccant ae 2354 0940 0124 2562 0956 2169 0956 
025-3 SCY Piserct ata oe o Ae — 0471 — 0173 — 0016 — 0544 — 0181 — 0431 — 0180 
dc, PASS MME Hote cs ayett one ene — 0562 — 0212 — 0019 — 0751 — 0277 — 0415 — 0117 
O%5-4 = AUP AG a8 otters oe oo eS — 0159 — 0072 — 0008 — 0324 — 0151 — 0059 0001 
8Cq Dies tere ape talas aie xeon vars agate — 1655 — 0722 — 0090 — 2181 — 0798 — 1046 — 0628 
284-5 SOL SIG creas elec yore oe -— 0550 —— 0241 — 0024 — 1113 — 0516 — 0180 — 0029 
8N;, = ON, — .03095 ......-.....- 1583 0762 0111 2168 0930 1097 0620 
dC, a WI ica OG 0 RTO EE SUT — 1085 — 0522 — 0082 — 1605 — 0608 — 0865 — 0513 
8cz = Fi, aed cae mocmOen par 0200 0120 0026 0261 0146 0137 0101 
Sly S— sie nabiow cure acd ur — 0034 — 0024 — 0007 — 0016 — 0010 — 0046 — 0036 
8(OPD) by differential Correction ....| — 0379 — 0144 0015 — 1543 — 0509 0361 0175 
Intermediate system. 
OPD* (by differential Correction)..... — 0251 — 0078 — 0033 — 0517 — 0163 — 0175 — 0002 
VED) SUCLELLACE) ciara siete othe iere svete io 3) <= — 0001 0011 — 0028 0343 0155 — 0323 — 0126 
82-1 MOC leet: ov stele) festa laa enn 0320 0150 0024 0324 0158 0324 0155 
O%5-3 Wb as Pehodic oa Oreck — 0383 — 0188 — 0028 — 0393 — 0212 — 0395 — 0201 
824-5 = bP Aninic. 5 ciel Os Opa iC — 0051 — 0024 — 0002 — 0112 — 0066 — 0003 0010 
dCg aa Mo ou tales oon ofSo Sa eres — 0212 —~ 0114 — 0020 — 0328 — 0159 — 0133 — 0082 
sNg= SN, — .0199........0.20-5. 0394 0242 0053 0502 0276 0324 0218 
8(OPD) by differential Correction.... 0068 0066 0027 — 0007 — 0003 0117 0100 
Modified System. 
OPD* (by differential Correction).... 0067 0077 — 0001 0336 0152 — 0206 — 0026 
OPD* (retraced). ....-. 6000s ereesees 0020 0050 — 0006 0347 0195 — 0309 — 0099 
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1.0 8L7--16.9ins $Z4.072 ins 


=-0074 ins 


8L’7=-2.94 8L=.719 
AC MARGINAL 85 NAN] Taxis 
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8LZ7=-9.05 S§L=.059 


-142 


\ 


Original 


system 


at the submarginal T-ray foci are nearly equal for 
H>0Oand H <0. For the marginal T-rays there is 
now a comatic under-correction, viz. O. P. D. more 
positive for the H< 0 than H>0. There is thus a 
zonal variation of coma which precludes its complete 
elimination. To the tangential focus (which includes 
both PerzvaL curvature and astigmatism dependant 
on the square of the aperture) again corresponds a 
straight line as reference sphere. This will have an 
inclination with respect to the y-axis dependant on 
the field curvature (obviously positive) of the objective. 
The asphericity of the extra-axial wavefronts is judged 
by their departure from such a line. Thus we note 
(fig. 8 a) as well as the residual coma, an overcorrec- 
tion dependant on even powers of the aperture. This 
however is tolerated since it assists the flattening 
of the effective field of the objective, while the worst 
of the comatic under correction at the margin, is 
eliminated by vignetting at the objective front. 

SL! (3L,) values (all measured from the Gaussian 
image plane) are quoted for comparison with the 
conventional method of characterising the correc- 
tion of an objective. 

In Table II (fig. 9) the technique of differential 
correction is applied to a4 mm. 85 N. A. dry achroma- 
tic objective. Such an objective is generally conside- 
red over-apertured and presents considerable difficul- 


Ne 
state 


2 7 
OPDy Wavelength) 


Modified 


system 


H<O 


ties to the designer (see Conrapy unpublished notes 
on microscope objective design). The objective chosen 
as the original system was severely under corrected 
at the zone, being corrected to a best focus considera- 
bly shorter than the paraxial focus. However, except 
for some over corrected coma, its field correction, 
notably field curvature, was not excessive. The main 
aim of the differential changes was thus to reduce the 
zonal under correction while maintaining as large a 
working clearance as possible. To be able to do this, 
the first requirement was the removal of the marginal 
over correction on axis. This was achieved in the first 
ten differential changes leading to the intermediate 
state (fig. 9 a). A large number of compensating 
changes had to be taken to maintain the convergence 
of the axial ray. Once a state of balance was reached 
at the edge of the aperture, leaving it in fact slightly 
undercorrected, it was possible to over correct further 
over the whole aperture and so remove to a large 
extent the zonal undercorrection. The last five chan- 
ges, calculated on the intermediate state, give the 
required result. The coma in the final form remains 
overcorrected but at only one third of its origina] 
value. 
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Appendix I. — An alternative method of deriving 
the focal shift terms used in the paper is given by 
the following construction. In figure 10a we consider 
the image (or object) space of figure 2, § III, alone 
and draw the bisector of the angle U. The perpendicu- 
lar to the axis at O* intersects it at M. O*M = o. The 
aberrations of an optical system are invariant not 
only at the ray focus, but also at any point on the 
bisector of the convergence angle (Hopkins, 1952). 
Hence the tranverse focal shift (normal to the axis), 
O*M, is equivalent to the longitudinal shift (along 
the axis) O* to O. From the geometry of the figure 

(S—t) u, U Uy sin U 
Sin. aad 2 Ag Deteircas UT win oe 
Substituting this in the transverse focal shift formula 
(BLASCHKE, 1954) we have 
sin U i, pe pine 
Loy Ge ly. |) 
A Cosabl Sf \ DA encos UE 
= — M(S — t) ® (sin VU, ut). 
This is the expression we have derived for the path 
N(B*B) of figure 2, § III. But in § IV, figure 3 the 
corresponding intercept is R* R, with reference spheres 
intersecting at P* instead of 0. This introduces a 
factor into the squared term to give 
sin U {sink U7 S— ty 
fee AE Ua i; 
4-++cos U 2 1+cosU S 
Nevertheless the error committed in § IV by setting 
R*R = B*B is quite negligible. 
We extend the treatment to the extra-axial case. 
In figure 10b PO, PO, -P*O*, P*O7, are respec- 


= NGS = 1) 


—-N(S—t) 


DIFFERENTIAL CORRECTION 
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tively pupil and T-rays through the original and chan- 
ged systems. 

PO; = S, POr =L, Por=t, Po, =t; 
< P*o. Pau, =. P*0, Pee 
the remaining notation being as before. We require 
the length of the normal to P*O* at 0; taken to its 
intersection M with the bisector of the angle PO,P(V). 
If Or M, is drawn parallel to P*O* then 


A ee V 
o = OM=(L—t) it —M,Or-+tan } Ci (u-+\ 


2 
Substituting 
eer (Da. 

M,0,=(L—t) u, tanV*+ | (S—t)u3—— aaa es 
and noting that 7 

7 Ve = V 

U* (u Ske 
we have 
o sin y(t, | sin V*-+cos V* tan (a, — x)\- 


es V 
— (S —?) uw, tan (a, a =) 


With the same approxiation as before 


V \ a" sin V 
2 4+ceosV * 


2 iu, 


4+cosV * 


tan (a, ast 


Hence 
Bieler on 
1+cos V aa 1-+-cos V | 
sin V 


+ (S —t) . 
( PANG + cos V 1+ cos V 


The focal shift term is given as in the axial] case by 
substituting c sin V* in the transverse focal shift 
formula. Collecting terms and making some further 
approximations consistent with those made in the 


osin V*¥=(L—t) u, ( 
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above, the result of general validity (but for a factor 
in the squared terms depending on the location of 
the reference spheres) becomes : 


N [\s 4) @ (sin V, u,) + 


Uy 


4-+cos V 
(1 — cos V) |. 


+(L —?) 


) Nik een ae ene 


® (sin V, u,) — 


Appendix II T-rays. — The formulae of the ray trace 
scheme for T-rays (extra-axial rays relative to pupil 
rays) are given in detail. The axial case follows imme- 
diately by setting I — 0, V =u —u = U,i.e. 
letting the pupil ray coincide with the axis. Similarly 
by setting sin J = 1,sin V = u,, 1. e. letting the aper- 
ture of the T-ray relative to the pupil ray become small 
the T-ray degenerates to a tangential fan. 

In figure 11a, C is the centre of curvature of the 
refracting surface 


<PCA=a, <ACA=«;AN=h,, Pg=Y,; 
AO =L, -AOy: = Ly POneS: 


With O; as centre a reference sphere is drawn through 
A intersecting PO in x. zy is the normal to AO, PJ is 


the parallel to AO, xP = «, < xPJ = V. 
Then PUTT eTeCe Tee rer eee eee rer er irre 
(i) Pascyieaeies sin V == h, a) fesin Ve 


CR is the normal from the centre of curvature a PO 
and produced meets AQ, drawn parallel to PO, in Q. 
Using r for the radius of curvature we have 


QR=rsinJ+rsin(U—«) =rsinI +rsin(V—J) 
or 


(11) r (sin J + sin V cos J —cos V sin J) = he 


A pouenene expression for ¢ is derived as rolloeae 

As in figure 2 § III we draw the reference sphere 
through P intersecting, in our more general case, 
the pupil ray at a distance « from A. The refracting 
surface itself is now considered a wavefront, viz. 
reference sphere, with axial focus C. We draw the line 
of invariant foci for such a wavefront i.e. the bisector 
of < PCA = « — a. O,, its intersection with AO, is 
the equivalent focus of the refracting surface on the 
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pupil ray. — «is then given by the focal shift QO, to O, 
along the pupil ray. (The negative sign is in accordance 
with the convention § III). 

_ In deriving an expression for <¢ no approximation is 
justified and we must return to (1) of § III. The place 
of u, is taken by 


see O.PO, => 


= | SE aes ere 


when —e = 8 (1 — ; 
sn(V+JI+4+4/) , 


sin V + sin oe) 


— e we have 


Replacing S by — 


sin 


s=— 


h., sin(V +1+7) 7 
sin V + sin (J + J) 


sin V 


5 2h, sin + 
(14 + cos V) cos (2 #7) + sin V sin u as !) 


For convenience of computation we multiply nume- 


Oe 


) and arrive at 


rator and denominator by 2 cos 


the expression 


amos amare Fane 4 ieee Re eae: 


e ,G=(1-+eosV)(cos I-+-e087)-+: 


+ sin V(sin J + sin J) 


Before considering the relations for the refraction of 
the ray we derive an alternative expression for Yv. 
It is used as a check on computations. Referring to 
figure 11 b, 
Yr = [Y— Y + r(cos «—cos«) tan U]cos u 
=r} sin(« — U) —sin(« — U) 
—r|sin(I + V)—sinI { 
Or in form convenient for computations. 
H PTEPUTTTETETETTPPTTT PEPE eeEEE LUTE e eee eee Cee etree Cece errr ere reer eee oe 
(iv) /E =r(sinI cos V+ cosJsin V),E=Yr +r sin L:. 


Two changes have to be considered in crossing the 
refracting surface Y,, to Y, and sin V to sin V’. (In 
the axial case only the latter). Referring to figure 411 b) 
Y,= [Pg — Ag tan J] cos I n > 
= r[cos I sin (x —«) — (1 — cos (« —«)) sin J]. 
Whence ” y 
N[Yr —rcosIsin(«—«a)] = 
= N'[Y!}—rcos I’ sin(«—«)] = 


= N sin I (cos (« — «) — 1), 
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the righthand side being invariant (invariant in the 
convenvional sense of a quantity which is the same 
a.r.and b. r.) Replacing mm this « —« by V + (J— J) 
we have 


(v) NY, = NY 
ir | sin V cos (I—T)-+ cos V sin (I—T) } A(N cos ) 


To derive sin V’ from sin V we use the invariance of 
the quantity 
N sin I cos (« —«) = 
= N sin J (cos V cos (I — J) —sin V sin(I—J)). 
To this we add and subtract 
N cos I sin V cos (I — J) and have the invariance of 


N cos (I — J) sin (I + V) — N cos J sin V. 


,K =N cos (I — Ne! 


(i) to (iv) with addition of the obvious results (4 and 2, 
denote consecutive surfaces) 
(vil) 


(viii) h,, = h,, —D sin Vi 


In the degenerate case of an axial ray they reduce 
to 


(ia) h= Y-+esia U 
(iia) hk =r (sin J + sin U) 
in I 
(iia)) e== 7 , G=(1+ cos UV) (4+ cos J) + 
+ sin U sin I 
N sin U—YK+YK’ 
(via) sin U’ = ww Ni a ,K=Neoos le 
(viia) (ees re ee 
(viiia) h, =h,—dsin U 
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The intensity at the surface of a semi-infinite 
diffuser with an infinite line source parallel to the surface 
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Summary. — Equations are developed for the intensity distribution at the surface of a semi-infinite isotropically scattering diffusing 
medium whose matrix may be of refractive index unity or greater, the source being a line lying on the surface ; numerical solu- 
tions are given for several values of the albedo for single scattering and hence of the reflectance. 

An extension of the argument leads to formal solutions for the intensity distribution for line sources below the surface. 


SomMAIRE. — On établit des relations qui donnent la distribution des intensités a la surface d@’un milieu semi-infini doué dune diffu- 
sion isotrope représentable par un modéle ott V indice est supérieur ou égal a Vunité ; la source est une droite de la surface. Des 
solutions numériques sont fournies pour plusieurs valeurs de Valbedo (correspondant a une diffusion unique) et par conséquent 


de la réflectance. 


Une généralisation du raisonnement conduit a des expressions de la distribution des intensités dans le cas ou la source linéaire 


est au-dessous de la surface. 


ZUSAMMENFASSUNG. — Es werden die Gleichungen fiir die Intensitdtsverteilung an der Oberfldche eines isotrop streuenden Mittels 
abgeleitet, das den unendlichen Halbraum erfiillt und dessen Matrix den Brechwert1 oder n hat. Die gradlinige Lichtwelle liegt auf 
der Oberjldche des Mediums. Die Zahlenergebnisse gelten fitr verschiedene Werle der Albedo des etnzelnen Streuvorganges und 
damit auch der Riickstrahlung. Eine Ausdehnung der Untersuchung fithrt zu den Formelwerten fiir die Intensitdtsverteilung bet 


einer linienférmigen Lichtquelle unter der Oberfldche. 


I. Introduction. — While many problems involving 
radiative transfer have been solved for plane parallel 
media and for media possessing spherical or cylindrical 
symmetry, such discussions have been restricted to 
cases involving intensity variations in only one dimen- 
sion. No solutions, however, have been given of pro- 
blems such as arise when a finite source is located near 
the surface of a diffuser, or when a diffuser is illumi- 
nated non-uniformly. 


In this paper we develop an approximate integral 
equation for the intensity distribution at the surface 
of a semi-infinite isotropically scattering diffuser for 
a line source on the surface ; and from the solution of 
this equation give formal solutions of the surface 
intensity for line sources parallel to the surface but 
within the medium. The method depends on a know- 
ledge of the intensity distribution for a line source in 
an infinite diffuser, recently derived by Jefferies [14] 
(1955). 

The present paper contains the relevant analysis 
and numerical solutions for semi-infinite diffusers 
having various values of the albedo for single scatte- 
ring. 


II. The intensity in an infinite isotropically scattering 
medium containing an infinite line source. — The 
model assumed for the diffusing medium is a 
matrix (whose refractive index V may be greater than 
unity) containing randomly spaced centres which 
scatter radiation isotropically ; radiation may also be 
absorbed by the medium. The coefficients of scatte- 
ring and “‘ true” absorption are denoted by o and « 
respectively, the extinction or attenuation coefficient 


otabyx;6=1—A= — is the albedo for single 


scattering. 

Let an infinite line source, emitting flux 4 m Q per 
unit length, be situated in an infinite diffusing me- 
dium. The direct flux from the source is converted 
into diffused flux by scattering, thus providing a 
source of diffused flux distributed through the medium. 
The source function B, which is the ratio of the diffu- 
sed flux produced in unit volume and emitted into 
unid solid angle to the extinction coefficient, is given 


by 


Q &o fe 
Deel == 
(2.1) yal 


e *% 


ea 

LV 2? — 6? 

JEFFERIES (1955) has shown that with an approxi- 
mation analogous to that of Epp1neTon for a plane 
parallel medium, the equation of radiative transfer 
may be written in cylindrical co-ordinates 


d2y ay 
2.2 | = 3 ¥(r1J—4 
(2.2) dott eee 3 ¥?(AJ —4 2 B) 


where ¢ is the perpendicular distance from the source 
and J the total intensity of diffused radiation, defined 
as-J = il I dQ where the integration is carried out 
over all solid angles and J dQ is the radiant power 
crossing unit projected area in an element of solid 
angle. JEFFERIES has given the solution of (2.2) in 
the form of tables of 
(J 9) 


6 xO 
for selected values of 6,. 
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Ill. The equation for the intensity at the surface of 
a semi-infinite diffuser having a matrix of refractive 
index unity and an infinite line source on the surface. 
Consider an infinite line source on the surface of a 
semi-infinite isotropically scattering diffuser emitting 
in all directions flux 4 Q per unit length, and let 
the total intensity of diffused light at the surface 
be Q. 9(r), where r is the distance from the source ; 
at any point on the surface the diffused radiation is, 
of course, confined to outward directions. Then if 
two such semi-infinite media be placed in contact 
enclosing the line source, the intensity distribution 
will become that in an infinite medium (§ II), denoted 
by Q.F(r). However, the latter may now be regarded 
as composed of (i) two components each independent 
of the presence of the other semi-infinite medium, 
totalling 2 Q¢(r), and (ii) components due to the dis- 
tributed or surface sources represented by the adja- 
cent surface of the other medium. 

Since the intensity of diffused radiation J is almost 
isotropic in an infinite medium, the flux passing from 
one of the above semi-infinite media to the other is 


Af ; ; 
closely GT per unit area, so that there is effectively a 
surface source over the interface emitting a total flux 
a ; 
of “> per unit area. The surface source gives rise to a 


different polar distribution of ‘‘ direct ” light from 
that of the line source but the effect of this difference 
is disregarded for the present. The “ direct ” light 
itself contributes nothing to the total intensity along 
the plane separating the two media, for there the sur- 
face source subtends zero angle. Thus at distance R 
from the source we find the approximate equation 


QM) = 20 eR) + | 
or 


1 1? 
Bl) 9(R)=> F(R) 5 / Fo) 6) dr. 


—o 


Since it is clear that o(R) = o(— R), the above equa- 
tion may be rewritten, with a slight change of variable, 


1 1 ss 
(3.2) (Rk) = 5 F(A) — go in F(R +r) 9(r) dr. 


IV. The intensity in an infinite medium with an 
infinite opaque strip source of narrow width. — The 
equations in § III have been derived neglecting ditfe- 
rences in the intensity distribution of diffused radia- 
tion originating from line and from opaque strip 
sources. For a source of non-uniform polar distribu- 
tion in an infinite medium the total intensity of diffu- 
sed radiation may be found by replacing the source 
with an equivalent source distributed throughout the 
volume and representing once-scattered radiation. 
The total intensity of diffused radiation at a point is 
thus given by the sum of the total intensities of the 


INTENSITY A THE SURFACE OF A SEMI-INFINITE DIFFUSER 25 


diffused and of the “‘ direct ”’ radiation from the dis- 
tributed source integrated throughout space. 

For an opaque strip source of narrow width emitting 
flux Q per unit length into an infinite medium it is 
not difficult to show in this way that at distance r 
from the source 


Ree ny ie 
Ai Sea I, 


Tt fmol 


ane | Fw rn 
2 Ki,(x 
4 SAGO | au ae 


where t? = r? + 62?— 2rov for a point in a plane 


containing the strip, or #2 =r? + p?—2roe Vi—v 
for a point in a plane at right angles to and containing 
the axis of the strip, and 
te ao hws cosh 
Kt (2) = | du . 
/ 0 cosh" w 


e 


The integrand in (4.1) is very small when xp > 1; 
thus the integral over e need not be extended beyond 
some value pe, > x". If r > e,, the values of ¢ and 
hence of J(r) are the same in both planes. It follows 
that at sufficient distances the intensity is independent 
of the polar distribution of direct radiation from the 
source, and this will probably be so at distances for 
which JEFFERIES’ solution of 


Fe) 
6m *Q 


approaches its asymptotic value 
e Fi care Cee ter 
2 108 G —1 ) He erry 
where @ = (3 a)”, 

At distances from the source small compared with 
x ', the total intensity of diffused radiation is also 
clearly independent of the polar distribution of direct 
radiation from the source. By inference, the variations 
at intermediate distances must also be small, parti- 
cularly for the larger values of 8. Reference to Ave 
fi 
the differences for different types of source are pro- 
bably not more than about 1 per cent, and are even 
smaller for 6 = 10; for 8 = 2 they may be appre- 
ciably greater. 

In view of the labour involved in evaluating (4.1), 
and the small differences expected, the results are 
assumed to be the same as for a line source for the 
purposes of the present discussion. 


FERIES’ solution suggests that for 8 = 5 1e., ) = 


\V. The equation for the intensity at the surface of a 
semi-infinite diffuser having a matrix of refractive index 
N and a line source on the surface. — Consider a semi- 
infinite diffuser with a line source lying in optical con- 
tact with the surface, and emitting into the medium 
flux 2 x Q per unit length. At distance r from the 
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source, and just within the surface, let the contribu- 
tion to the total One from radiation directed 
towards the surface be Q 9(r) ; thus the flux incident 


(7) 


per unit area. Since the transmission of the surface 
to uniformly diffused radiation incident from within 
is J/N2, where J is the transmission of the surface to 
uniformly diffused radiation incident from air, the 


on the surface from below is very closely a 


escaping flux is very closely e Q o(r) (5/N?) per unit 
area. 

Now let two such semi-infinite media be placed 
almost in contact, the source emitting flux 2 x Q per 
unit length into each medium. The distribution QF(r) 
is thus restored. As before, this may be regarded as 
composed of (i) components independent of the pre- 
sence of the other semi-infinite medium, and (ii) com- 
ponents due to the distributed or surface source repre- 
sented by the adjacent surface of the other medium. 
In the case of (i), radiation reflected at the surface 
and transmitted through it are complementary, so 
that the total intensity due to these components is 
2 Q 9(r). The surface source due to radiation transmit- 
ted from one medium to the other effectively emits a 


total flux of 5 Q(5/N?) F 


ding the ere of ate reflections between the 
surfaces. Any element of this effective plane source 
has a polar distribution of radiation which depends 
on the glossiness or optical state of the surface but, 
as before, it is sufficient to assume that the total inten- 
sity of diffused radiation due to a source of this type 
is the same as for a distribution of line sources. Thus 
at distance R, 
aif 


Q FR) =2@ oR) + 7 


or 
(5-4) o(R) = > F(R) atl 


) per unit area, this inclu- 


r) oR — r) dr 


F(R-+r) 9(r) dr. 


This equation reduces to (3.2) when V = 1. 


VI. Integral properties of o(r) and F(r). — Since 
when NV = 1 the escaping flux is very closely eA per 


unit area, it follows that for diffused flux 2 x @Q incident 
on the medium the total reflected flux is 


in g(r) dr. 


The ratio of the latter to the former is the reflec- 
tance for diffused incident radiation, R so that 


(6.4) [ o(r) dr=40R. 


[Orr. Acta 


Integrating (3.1) with respect to R over the range 
—o too 


«foe 9(R) dR = » fa F(R) dR — 
oe ae R—r)drdk. 


Carrying out the double integration over R first, it 
follows from (6.1) that 


a Sr R 


When WN > 4, it may be shown that 
sa 2 
[ g(r) dr = ae R, 


where ‘A is now the reflectance for diffused incident 
radiation omitting FREsNEL reflectance, or 


where ®p is the total reflectance for diffused incident 
radiation. 


(6.2) 


(6.3 


R= Rp cao 


VII. Solution of the intensity equation. — Equa- 
tion (3.2) may be recast slightly by noting that 9(r) 
may be expected to decrease with increasing r more 
rapidly than does F(r). Thus 


(7) 


xr) = Fr) 


is expected to vary less rapidly than 9(r). Equation 
(3.2) may then be rewritten as an equation in y(r) : 


1 


(7.4) (RB) =5— fF (+1) F(r) x(p) dr 


Ew 


The function y(r) is now replaced by a polynominal 
in 7; or alternatively y(r) is allocated unknown 
values «, @, y, 8... at equally spaced intervals of r, 
intermediate values being obtained in terms of these 
by Lagrangian interpolation. For a given value of wp, 
F(R) is obtained from JerreriEs’ tables and by no- 
ting that for large values of R, 


lined rye 
Baar gee Goan 


where 8 = (3 4) 4/?; graphical interpolation for small 
values of R enables values of F(R) to be obtained 
with sufficient accuracy at close intervals. The inte- 
gral in (7.4) can thus be Rake numerically as a 
linear function of «, 8, y, 8,..., for given values of R; 
for the latter, a choice of ate same values as for Me: 
6, etc., enables (7.1) to be replaced by a set of linear 
simultaneous equations for «, 8, y, 8,..., which can be 
solved in the usual way. 


(32 xR) , 
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In this method of solution, two points need special 
mention ; (i) it is desirable to replace the variable r 


in the integral of (7.1) by y = \/3 Axr, in terms of 
which F(R) is more readily expressed ; (ii) rather than 
evaluate the integral for values of r for which ,(r) 
would need to be extrapolated, the contribution for 
these larger values of r (which will be small in any 
case) is found by assuming an exponential decrease 
in the integrand, whose rate is estimated from the 
ratio 8/y obtained by trial solutions. 


The solutions given in Tables 1-3 have been obtai- 
ned with four values of y(r), the corresponding values 


of y and we being indicated in the appropriate table, 
while graphs of 2”) in figure 1 represent the functions 
x 
derived using four-point Lagrangian interpolation 
between the y(r). 
At no stage can the accuracy of these results be 
expected to be better than one per cent, and it may be 


several times poorer. The EppINGTon approximation 


y (@,=299 / 300) 
6-0 01.02 03 O4 O05 086 O07 08 
‘ Pad | 


5:0 


4-0 


O——O @,=299/300 
x—*x @®,=74/75 
A—A @®,=11/12 


2:0 


1:0 


ime G02 e000 20-6) 906) 410 | 12m ident 
y (@,= 74/75 & 11/12) 


or 
Fic. 1. — The function ae The values of the abs- 


es 299 
cissae, y = V 3A xr, depend on ég, those for 6) = >— 


74 
being indicated at the top, and for a» = 7 and DB 
at the bottom, of the graph. Light and heavy curves 
are for N = 1 and N = 1.5 respectively. 


TABLE 1 


Intensity distribution at the surface of a semi-infinite isotropic 
diffuser with a line source on the surface 


299 
9 = 559 = 0:996667 
¥y 
0 0.25 0.5 0.75 
F 
ra ©) ol485 15419 0927 ee O61s 
6 ay % 
( A(r) 0.2074 0.0573 0.0323 0.0176 
Med a o(r) 3.904 0.528 0.179 0.064 
(234 
( x(r) 0.2865 0.1493 0.1012 0.0728 
Newaeo 1 _ 
(na o(r) 5.393 1.376 0.561 0.267 
foe 5 


TABLE 2 


Intensity distribution at the surface of a semi-infinite isotropic 
diffuser with a line source on the surface 


74 
®o = zz = 0.986667 


2.554 0.885 0.427 0.217 


(207) 0.2623 0.0879 0.0510 0.0266 
Nd Sie LTE 3.966 0.460 0.129 0.034 
ie 

\ 
(40 0.3462 0.2029 0.1485 0.1067 
/ 
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TABLE 3 


Intensity distribution at the surface of a semi-infinite isotropic 
diffuser with a line source on the surface 


11 
y= Fy = 0.916667 
y 
0 0.5 4.0 a5 
zt 1.792 0.843 0.430 0.234 
6 &o% 
- X(r) 0.3533 0.2494 0.1863 0.1462 
ak ja o(r) | 3.482 1.156 0.441 0.188 
x 


(tr) 0.4207 0.3598 0.3160 0.2816 
5) 4146 1.669 0.747 0.362 
x 


is known to introduce errors of this magnitude, and 
JEFFERIES’ equation of transfer probably introduces 
comparable errors. Further a small error has been 
introduced in deriving (3.1) by assuming that the 
line source and the surface source produce intensity 
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distributions along the plane of separation of the 
two media which are of identical form 9(r). Finally, 
the solution depends on the replacement of y(r) by a 
polynomial, and this necessarily introduces further 
errors, the effect of which may be noted in the compu- 


¢(r) 


ted value of &~ at the source, which is nearly 2 per 
% 74 : 299 
cent larger for é = 75 than for 6) = 300? whereas 
it should be smaller. However, extremely useful 
checks can be made via the integral relations (6.1) 
ake) 
Fir) dr obtai- 


—_—oA 


and (6.2). The numerical value of | 


ned from JrerrenrEs’ solution agrees to within one per 
nih 
1— A 
rences in polar distribution of the sources) for the two 
higher values of 6 , but is rather higher for the 


ran 


cent with the value (derived neglecting the diffe- 


lowest. value of 6». Further the integral | o(r) dr 
evaluated numerically from our solutions agrees satis- 
factorily with the value 4‘, as may be seen from 


Table 4. 


TABLE 4 
ice) (oe) 
The Integrals F(r) dr and iii o(r) dr; N=1 
x 0 eee 
; 11 74 299 
,° 12 75 300 
RK IB GYOUS (HAG Shoo Sewanee 0.5111] 0.7668) 0.8756 
“* ? On EppineaTon approx. | 0.5199] 0.7695] 0.8764 
toa) 
jel(a Wales nes ae: 14.11 |41.7 |89.6 
0 
AW WCAC Ee, seca t, 13.14 |41.3  |88.4 
1—R ) On Eppineton approx. |13.61 |42.0 89.1 
oo 
o(r) dr 6.65 9.60 {10.77 
—O 
Eexacts (ee is ae sistas soa ecowe 6.42 9.64 |11.00 
4nf— 
On EDDINGTON approx... | 6.53 SLO 7a AC On 


VIII. The intensity distribution at the surface of 
a semi-infinite diffuser with an infinite line source 
beneath the surface. — For an infinite line source at 
distance Z beneath the surface of a semi-infinite 
diffuser, and emitting flux 4 Q per unit length, let 
the total intensity of diffused radiation at the surface 
be Q. 6(r) at distance r from the projection of the 
source on the surface. 

For a matrix of refractive index unity, suppose as 
before that a similar semi-infinite diffuser is now pla- 
ced in contact with the surface, the intensity distri- 
bution of the diffused radiation being thus increased 
to QF(¢), where ¢ is the distance from the line source. 
At the interface this total intensity is composed of 


[ Orr. Acta 


(i) the independent component due to the first me- 
dium (ii) the component in each medium due to the 
distributed or surface source represented by the adja- 
cent surface of the other medium, and (ili) the com- 
ponent in the second medium due to irradiation by 
residual direct light from the source, which would nor- 
mally have escaped from the first medium. The last 
component is somewhat different in character from 
those discussed previously, and will now be evaluated. 


FIG. 2 


Fic. 2. — Line source AB beneath the surface of a semi-infinite 
diffuser C D L M. 


Let AB (fig. 2) be the line source at distance Z below 
the surface CDLM of the diffuser. Then the surface 
source due to direct radiation on the upper medium 
can be evaluated by finding the direct illumination 
on the surface CDLM. From an element of source of 
length dl at G, the illumination d# at H, distant r 
from the projection of the source on the surface, is 


OCC a 


2 


(ee COS L 

where x = GH is the distance of the element of source 
from the point under consideration, and 7 is the angle 
GH makes with the normal to the surface. Thus 


co 53 3 
6 =a Qe cost 4, 
—n ie 
f Z 
But cos 7 = A and Id! = «dz, so 
ae e%* 
(8.1). Bt) == 21052 dx. 


“6 2/x2— Z2—r2 
This is the flux entering unit area of the upper medium: 
assuming here too that the surface source is associa- 
ted with a 9(r) of the same form as the line source, at 
distance R from the projection of the source on the 
surface it gives rise in the upper medium to a compo- 
nent of the total intensity of diffused radiation 

4 (ee) 

Dieses 

We may thus proceed to write down the approxi- 

mate equation for 6(R) : 


E(r) o(R —r) dr. 


na 
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Q F/R? + Z2) = 
R) + ee / _ F(Vr?+22) o(R—n) dr + 


1 (o/¢) 
a ae is E(r) o(R —r) dr. 


The formal solution for 6(R) follows immediately : 


(8.3) = F(R? + Z*) — 
di (oe) 
-n/ F(Vr? + Z2) @(R—r) dr — 
i\ = Er) 
= 6 o(R —r) dr. 
For sources well below the surface, so that Curis 


negligible, E(r) is negligible and (8.3) simplifies to 
(8.4) OR) = 


F(\/R? + Z2) — = [- F(Vr? + Z?) «(R—r) dr. 


It must be remembered that (8.3) represents only 
the diffused radiation. The total radiation field is 
found by adding the direct component. 

For a matrix of refractive index JV, suppose that an 
identical semi-infinite medium containing a similar 
line source is placed almost in contact with the first 
medium, the line joining the sources being perpendi- 
cular to the surface. Then the total diffused intensity 
at any point on the surface distant R from the line 
joining the sources is 2 QO F(\/R? + 2). This is com- 
posed of 

(i) two independent components totalling 2 Q0(R) 
where Q6(R) is the total outward intensity in a hemis- 
phere just below the surface when only one semi-infi- 
nite medium is present ; 

(ii) two additional components totalling 


oe de E(r) t(r) o(R —r) dr, 


2 
arising from direct radiation from the one source inci- 
dent on the other medium, t(r) being the transmit- 
tance of a single surface of refractive index J for 
‘“‘ direct ” light from this source — in the case of a 
polished surface, t(r) being given by FResnev’s well- 
known rules ; and 

(iii) components due to the distributed or surface 
source represented by the adjacent surface of the 
other medium, and totalling 


DODe pe —— 
om | PWR +2) te 
(see § IV). 
Thus 


20 F(\/R?+ Z) = 


5 a ——— 
~=200(R) + ae P P(r +24) o(R—2) dr + 


+o fo 


The formal solution for 6(A) follows : 


) 0(R —r) dr. 
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= F(\/ Re + Z2) — 
F (Vr? + Z*) o(R—r) dr — 


(8.5) 


‘ Ds 

ae onde 
1 ieee 

= i . ss e(R —r) dr. 


For sources well below the surface, so that e “7 is 
negligible, E(r) is negligible and (8.5) simplifies to 


(8.6) (R) = F/R? + Z2) — 
J 2 
ae aft F(\/r? + Z?) o(R—r) dr. 


Again, 0(R) is only the Gutneed component of the 
outward radiation field immediately below the sur- 


face ; the total diffused intensity is o(R)| 2 = <n) ; 


the total radiation field is found by adding the direct 
component. 


IX. Discussion. — It has been pointed out that the 
function ¢(r) depends on the polar distribution of 
direct radiation from the source. We have shown that 
as 6» approaches unity the solution becomes almost 
independent of the polar distribution, so that the 
results given here are good approximations for pro- 
blems such as the distribution of intensity along the 
surface of a semi-infinite diffuser for light incident 
from any cone of illumination. The resolving power 
of the surface used as a projection screen may be assess- 
ed from the range + R in which half the light is 
reflected. In Table 5 this is given in, terms of 


Y= 3axR. 


TABLE 5 


Range + Y in which half the light is reflected 


Y 
Og 
Viel N= 1.5 
ay 0.10 0.17 
| 300 . . 
| ie 0.18 0.26 
75 a — 
- 0.35 0.40 
42 5 


The above results also give, with appropriate 
changes of terminology, the distribution of thermal 
neutrons along the surface of a semi-infinite medium 
with a line source of thermal neutrons parallel to the 
surface ; for the diffusion equations governing radia- 
tion and thermal neutrons are formally identical. 
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Gekoppelte Interferenzfilter 


Urricn GRADMANN 
Institut fiir angewandte Physik, Marburg. 


ZUSAMMENFASSUNG. Einfache Interferenzfilter bestehen aus zwei parallelen, in Dielektrika eingebelteten Metallschichten. Geht 
man zu Systemen mit mehr als zwei parallelen Metallschichten iiber, setzt also mehrere solche Einfachfilter hintereinander, so 
multiplizieren sich deren Durchlassigkeiten nicht einfach, es treten vielmehr Kopplungseffekte auf. Fir Filter mit 2, 3 und 4 
Metallschichten werden Formeln fiir den spektralen Verlauf der Durchldssigkeiten angegeben und diskutiert. Numerisch gerech- 
tete Beispiele werden in Diagrammen dargestelilt. 


SomMArRE. — Les filtres interférentiels simples se composent de deux couches métalliques paralléles entourées par des diélectriques. 
Sil’on passe a des systémes comprenant plus de deux couches métalliques paralléles, c’est-d-dire lorsqu’on met en série plusieurs 
filtres simples leurs intensités ne se multinlient pas simplement ;,de plus, des effets de couplage apparaissent. On donne et 
discute les formules pour la variation spectrale des facteurs de transmission pour des filtres a 2, 3 ef 4 couches métalliques. Des 
exemples calculés numériquement sont représentés sur des diagrammes. 


Summary. — Simple interference filters consist of two parallel metallic layers alternating with dielectrics. For systems of more than 
two parallel metallic lavers, i. e. when several simple fillers are placed in series, their transmissions are not simply multiplied, 
but coupling effects appear. Formulae are given and discussed for the spectral variation of the transmission factors for filters 
of 2, 3 and 4 metallic layers. Diagrams are given of calculated numerical examples. 


1. Einleitung. — Reflexion und Durchlassigkeit und die durchgehende 
eines einfachen Interferenzfilters aus zwei diinnen i (+ 72° 
Metallschichten und einem dazwischen liegenden Die- D=De as c 


lektrikum sind theoretisch seit langem geklart. In 
seiner Patentschrift zu diesem Thema zeigt Grrr- 
cKEN [1], wie man durch Hintereinanderschalten 
entkoppelter, d. h. unter einem spitzen Winkel zuein- 
ander stehender Elementarfilter zu Filtersystemen 
mit sehr scharfen Linien kommen kann. Werden die 
Filter jedoch nicht entkoppelt, stehen sie also parallel, 
so treten kompliziertere Verhaltnisse auf, die hier 
durch Zuriickgehen auf die Maxwe.v’sche Theorie 


diskutiert werden sollen. Das Problem steht in Ana- Ass. 1, — Ein System paralleler Schichten, be- 
Jogie zur Vierpoltheorie der Hochfrequenzschaltungen, sherry een Drechinn Bateae aad sige 
worauf ScHUSTER [2] hingewiesen hat, und man erhalt bettet. Einfallende ((), reflektierte (9) und 
dementsprechend auch hier Bandfilter, d.h. Filter mit durchgehende Welle (2). 
einem verbreiterten Durchlassbereich hoher Flan- 
kensteilheit (1) [3]. Wir setzen dann 

Unter einigen vereinfachenden Voraussetzungen, D NR 
die das Wesentliche klar herausstellen, sollen nun (la) == p% “und (1b) == 4 
solche gekoppelten Filter besprochen und Zahlen- Go Eo 


beispiele angegeben werden. Wir beschranken uns und analog fiir senkrechten Einfall von unten, also 
auf senkrechten Einfall, nehmen alle Metallschichten ‘ ‘ K 
als gleich an, betten sie in ein einheitliches Dielek- mit der z— Richtung, — = bd” und = = 1, 


trikum ein und vernachlassigen die Dispersion der. Xo , 
optischen Konstanten. ferner fiir die Teilsysteme im gleichen Sinne 9°, 6° 
2. Grundlagen. — Ein System paralleler Schichten usw. Dabei soll die Null als oberer Index bedeuten, 


sei gegeben, bestehend aus den Teilsystemen I und II, ae one Beis aul — Sean desselben hoon 
als ganzes eingebettet in ein Dielektrikum vom Bre- Romans bai ake ct Stet ete lial) 


chungsindex ny. Wir betrachten nun eine senkrecht br Dd, 
einfallende, der z— Richtung entgegenlaufende, ebeen (2) bear {2050 
und linear polarisierte Welle (Abb. 1). oe 
way = tr tered) yar ee 
& = & e o(it ¢ ) ? (3) ” F ; a re : “ ) 
Die am Gesamtsystem reflektierte Welle sei et 
- 3. Das Einfachfilter. — Das in diesem Abschnitt 
io(™ behandelte Einfachfilter besteht aus zwei Metall- 


schichten der Dicke h, vom Brechungsindex n, die im 


(1) Inzwischen wurd Iche Filt <peri : F ape , 
Gemeiekneit [3]. en solche Filter auch experimentell (*) Siehe z.B. den Artikel « Optik dinner Schichten » in 


der Neuauflage des Handbuchs der Physik (im Druck). 
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Hy (= 1033 mp) 


b) 


Ze 


h(=30mp4 


H(= 884 m4) 


I H, (= 1033my4 


GEKOPPELTE INTERFERENZFILTER 


ABB. 2. — Verschiedene Filter. 


a) Einfachilter, 6) Doppelfilter, c) 

Dreifach filter. 

Erste Spalte : Aufbau der Filter. 

Zweite Spalte: Nenner der Durch- 
lassigkeiten,;f (a) bzw. oh (b) in 
der Zahlenebene. Die offene Kurve 
in (b) ist der erste Nennersummand 
von »p, der Kreis das Kopplungs- 
glied, die Pfeile sind der Nenner 
selbst. Parameter 9 bzw. 9, inGrad. 


Dritte Spalte : Durchgelassene 
Intensitaten in % der eingefallenen. 
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Abstand # in ein Dielektrikum vom Brechungsindex 
ny eingebettet sind. (Abb. 2a). 

Zuniachst seien Reflexion or (!) und Durchlassig- 
keit )d einer einzelnen Metallschicht berechnet, die 
im Sinne von Gl. (1) definiert und auf z = 0 bezogen 
seien. Man geht dazu mit den bekannten Fres- 
NEL’schen Koeffizienten y und ) in (2) und (3) ein (die 
‘« Systeme ’ I und II sind hier die beiden Grenzfla- 
chen), beachtet aber dabei, dass diese Koeffizienten 
auf die jeweiligen Grenzflichen selbst bezogen sind. 
Sie stehen zu unseren Gréssen also in den folgenden 
Beziehungen (Abb. 3 b) : 


TD. e012 

eee 2m — D(%r) eed Se — p09 elf1—Po)% 
To + ty G4) E, e°o7 

5 he Ilo ie (21) Mo erie = 7? e2 2071 
No tn,  €(%) E, eH 


b) ~ ZX 


A 


Ass. 3. — a) Eine Metallschicht der Dicke h 
vom Brechungsindex n,, in ein Dielektrikum 


vom Brechungsindex ny eingebettet. b) Grenzfla- 
che zweier Medien mit den Brechungsindizes 
Mo und ny. 


Damit erhadlt man 


een | 
Notun 
(4) b= oe eroh — jhe etl 
Ne — th, 
1— e- 2 Pm” 
Wo => 1, 
No —— Th 4 —e—2Pmh 
Sle es yee ihe 
v Be fee Se? ome 
Ny +1, 
wobei 
2 nin; 
(6) 9, = ‘und a die Vakuumwellenlange ist. 


Im Rechenbeispiel ist die Schicht als 30 mu dicke 


Silberschicht, ferner n) = 2 angenommen. Der Zah- 
lenwert fiir n,, nach Kraurkrimer, Murmann, WoL- 
TER wurde dem Buche von H. Mayer [4] tiber diinne 
Schichten entnommen. Fiir h = 30 mu, 4 = 600 mu 
findet man dort 


(4) Der Index wurde vorgestellt, um Konflikte mit den 
Bezeichnungen des oben zitierten Handbuchartikels zu ver- 
meiden. 
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Setzt man diesen Wert in (4) und (5) ein so erhalt 
man 


(4 a) od* = 0,405 e826! 
(5 a) ot = 0,888 @7)1801 


Dass sich 9dD* und or in der Phase um nahezu 90° 
unterscheiden, folgt aus (4) und (5) unter Beachtung 
onsen 
die Phase einer komplexen Zahl 2). 

ot und ,) zeigen an sich natiirlich Dispersion, 
wurden den weiteren Ausfiihrungen jedoch als Kon- 
stanten zugrunde gelegt, um das Grundsatzliche klarer 
herauszustellen. 

Um nun Reflexion ,r und Durchlassigkeit ,b des 
Einfachfilters zu bekommen, die im Sinne von Gl. (1) 
auf z = 0 bezogen seien (Abb. 2a), hat man nur mit 
(4) und (5) in (2) und (3) einzugehen, wobei freilich 
die Koordinatentransformation fiir die zweite Metall- 
schicht zu beachten ist. Man erhalt 


(od)? 
4 — por’ ot e2 0071 


a — = und wird weiter unten benutzt (¢, sei 


(od)? 


(Po (or)? e2 09(H—h) 


1 — e200! } (gx)? e? Pot — (gp)? | 


b= 


we ==" We 


aaa Be e2 0071 


A toes e-2 Q9(H—h) } (G 12 — 


(ob*)? { 


= of 


4 — (or)? e2 Q9(H—h) 


wenn man or’ = or e”°o" beachtet. Mit 


6) pees ee) Se = 
i ot 
und wenn man annimmt, dass sich »)* und or um 90° 
in der Phase unterscheiden, was, wie oben bemerkt, 
hier mit hinreichender Naherung zutrifft, wird aus 
den letzten Gleichungen 


(ob)? 
Z AD = - 
(7) Q 1—| or 2? e739 
| ety eae ot 2 od (2 
ay theta ef?) | oft | ob 7 | 


do eerrtieme 


Da od und or Konstanten sein sollen, stellt der 
Nenner von (7) und (8) einen Kreis dar, wie in Abb. 2a 
zu sehen ist, aus der man auch deutlich sieht, wie die 
Maxima fir ein Durchlassigkeitsfilter umso schirfer 
sind, je mehr sich »y dem Betrage nach 14 nihert. 
Daher wurde die Silberschicht 30 my stark gewahlt, 
was einigermassen scharfe Maxima ergab. 

Als Spiegelung des Nenners am Kreise | r |= | 9b |? 
stellt auch ,b einen Kreis dar. Dasselbe gilt auch 
fiir iv, wie man leicht nachrechnen kann. Abb. 4a und 
4b zeigen ,y und ,) in der Zahlenebene. 


iv und ,b sind in 27 periodische Funktionen von 


vol. 8, n°. 1, March 1956] 


05i 


-05i 


_ Ass. 4. — a) Reflexion jy und b) Durchlas- 
sigkeit ,) des Einfachfilters in der Zahlenebene. 


und damit auch periodische Funktionen der Fre- 
quenz v. 
Maxima der Durchlassigkeit treten auf fiir 


(9) 


” 
! 
\ 
bd 
é 
40 
20 
0 
-20 
qa. 4 
Ass. 5. — Amplitude und Phase von a) ,v und b) ine 
4 
c 
a 
& 
a +4 %e 
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Scorn 
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& 
Bo iP 
a 44 47 50 53 56:10""Hz 59 v—> 
Ass. 6. — Ausschnitte aus den Durchlassigkeits- 


spektren zweier Filter : 
a) Filter 1: k = 5; H, = 884 mu. 
Maxima von |,p| fir v= 5,00 ; 5,87 ; 6,74 usw. 


xelOt2 secs: 
Halbwertsbreite (bei 5 x 104 sec!) 8 mu. 
Zehntelwertsbreite 26 mu. 


b) Filter Il: k = 6; H,, = 1033 mu. 
Maxima von |,f| fir v = 5,00 ; 5,73 ; 6,46 usw. 


MOS tee 
Halbwertsbreite (bei 5 x 10! sect) 7 my.. 
Zehntelwertsbreite 22 my. 


k heisse “ Ordnung ” des Maximums. Ist eine 
auszufilternde Wellenlange gegeben, so dient (9) 
zur Bestimmung des Abstandes der Metallschichten. 
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Die Ordnung k ist dabei noch frei waihlbar und _ be- 
stimmt nach Gl. (9) die Breite der Linie, daja ,r und ,b 
Funktionen von @ allein sind. Erhoht man bei festem 
4 die Ordnung k, indem man den Abstand H der Me- 
tallschichten vergréssert, so erhdlt man _ scharfere 
Linien. Diese riicken dabei niher zusammen (sie sind 
in der v-Skala aquidistant, s. Abb. 6). 

Als Beispiele sind zwei Durchlassigkeitsfilter fiir 
4 = 600 my, v= 5 x 101*sec-!,k = 5 und k = 6 be- 
rechnet. Abb. 2a und 4-6 stellen die Ergebnisse der 
Rechnung graphisch dar. Abb. 2a, 4 und 5 entspre- 
chen numerisch dem Filter I von Abb. 6. 


4. Das Doppelfilter. — Wir betrachten ein Dop- 
pelfilter aus 3 identischen Metallschichten nach Ab- 
schnitt 3 in den Abstanden H, und H,, (Abb. 2b). Wir 
betrachten es als Durchlassigkeitsfilter, berechnen 
also nur seine Durchlassigkeit ,b (im Sinne von Gl. (1), 
bezogen auf z = 0). 

Um Formel (2) anwenden zu kénnen, fassen wir 
zwei Schichten zum System I zusammen und fassen 
die dritte als System II auf (Abb. 2b). Dann ist 


FO 
atin 
= 1d od 
4 — yx e2 P0lHi—*) oy @ 2 Pol + Ai) 
oe ib od 
1 — av ot e290) 
(ob)* 
10 = 
ie ° = WN, +R 
mit 
N,, = 1 —|orP ei Pk Kiley |2| od P e—ile,+9,,) 


[Bedeutung von 9, siehe Formel (6)]. 

Ohne K im Nenner und bis auf einen Faktor (bd)? 
stellt (10) offensichtlich das Produkt der Durchlis- 
sigkeiten von zwei Einfachfiltern im Sinne von 
Abschnitt 3 dar. K ist also eine Art Kopplungsglied. 1) 
In der Gauss’schen Zahlenebene ist K ein Kreis. 

Gibt man den beiden Elementarfiltern zusammen- 
fallende Maxima, so erhilt man bei der betreffenden 
Frequenz auch fir das Doppelfilter ein Hauptmaxi- 
mum, das jedoch durch das Kopplungsglied aus einer 
scharfen Linie zu einer Bandfilterkurve verformt 
wird, wie man sie ahnlich in der Hochfrequenztechnik 
kennt. Abb. 2b zeigt das in der Gauss’ schen Zahlen- 
ebene und stellt zugleich ein Rechenbeispiel gra- 
phisch dar (H; und Hy gleich denen von Abb. 6, 
Hauptmaximum bei 4 = 600 mu). Die Halbwerts- 
breite betragt 17 my, die Zehntelwertsbreite 24 mu. 
Man hat also bei einer gewissen Verbreiterung eine 
wesentlich hohere Flankensteilheit im Vergleich zum 
Einfachfilter. 

Storend mag in manchen Fallen das tiefe Minimum 


1) Das graphische Verfahren wurde von WOLTER in einer 
Vorlesung WS 1954/55 angegeben [5]. 


34 U. GRADMANN [Orr. acra 


in der Bandmitte sein. Um es zu beseitigen, muss man } 14 —| ot |? e719 bleibt dann unter 45°, die von 
die geschwungene Kurve in Abb. 2b, d. h. den ersten sale ewan 
Nennersummanden von Gl. (10), mehr in die positive 1i== Jing jaue unter : 


1Si 


7a) h= 30m pa 
(durchgel. Jntensitat 
d. Einzelschicht> 16%) 


To 


durchgel. Jntensitat gFi ilters 


z 
z | 
RE yp 
-120 -60 0 60° 120 
'H 
g 
t 
HO, 
& 
% 
0 Si S 
7b) h= 20 mp RS 
(durchgel. Jntensiat ~ 207 
d. Einzelschicht 30%) & 
o 
-120 -60 0 60° 120 if 
g 
it 
U 
8 
7) 
05i o 
7c) h = 15mpa ; Eady, 
(durchgel. Jntensitat i? 
a. Einzelshicht 40%) 3 
& 
p 
-120—— -60 0 60° 120 


Azss, 7. — Doppelfilter mit verschieden starker Kopplung. Erste Spalte : Aufbau des Filters. 
Zweite Spalte : Der Nenner von ,)* (s. Gl. (11)) in der Zahlenebene. 


Dritte Spalte : Durchgelassene Intensitaten in % der eingefallenen. 


Halbebene ziehen und die Eindellung bei » = 0 ein- Die Verhaltnisse sind in Abb. 7 dargestellt. 
ebnen. Das gelingt aber offensichtlich fiir kleinere | gr |, oe 
wie man aus Abb, 2a ablesen kann : Die Phase von Fir a) h = 30 my, ny, = 0,12 — 4,15 i, 
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Deb 20 me, 1, = 044— 4201. 
c) h = 15 my, n,, = 0,14 — 4,25 i. 
(Werte aus Mayer I. c.) ist. 


(od)* 


(11) : : 
feces | ot 2e—i9 {2 ais ov |? | 4d \2e—2i9 


b* ee 


nach gerechneten Werten als Funktion von © darge™ 
stellt. Zur Vereinfachung wurden also die beiden Ele- 
mentarfilter gleich gemacht (9, = 9,, = 9). 

Abb. 7 zeigt, wie die Kurven sich fiir diinnere Me- 
tallschichten ins Positive ziehen ; daneben gibt sie den 
Intensititsverlauf wider. 

Die Eindellung ist fiir kh = 20 my nur noch ganz 
gering. Die Einebnung wird durch Linienverbreite- 
rung und geringere Flankensteilheit erkauft. 


5. Das Dreifachfilter. — Endlich sei ein Dreifach- 
filter besprochen. Es besteht aus 4 Metallschichten 
im Sinne von Abschnitt 3 in den Abstanden Ay, Au: 
und Hi: (Abb. 2c). Wir fragen wieder nur nach der 
Durchlassigkeit ,d, gehen dazu auf Gl..(2) zuriick und 
fassen je zwei Metallschichten zu den Systemen I 
und II zusammen, deren Durchlassigkeit Gl. (7) angibt. 
Dann ist die Gesamtdurchlassigkeit 


Eine elementare Zwischenrechnung ergibt 


tam oO 


oa HE ae Nee ae 
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mit 
Ni, aoa 1 — Vetsle ein | ; i — | ov |? | od [2 ein 


} eter 4 ete — [| ob [2+ 2] or|2] elerten! 


Auch hier kann man in K* das Kopplungsglied 
erkennen. Der ganze Nenner ist aber so kompliziert 
gebaut, dass auf seine Darstellung in der Zahlenebene 
verzichtet sei. Zur numerischen Berechnung eines 
Beispiels geht man besser von (12) aus. Auf diesem 
Wege ergab sich mit Hin = Hi, Hu von Abb. 6, 
die Kurve der Abb. 2c. 

Auch fiir das Dreifachfilter ergibt sich, wie zu er- 
warten, eine Bandfilterkurve mit hoher Flanken- 
steilheit, das Band erscheint jedoch gegeniiber dem 
Doppelfilter schon etwas geglittet. Im sichtbaren 
Teil des Spektrums ergab eine graphische Abschatzung 
Nebenmaxima mit bis zu 17% der Hauptmaximums- 
intensitat. 

Durch Variation von Dicke und Abstand der Metall- 
schichten, durch Wahl anderer Metalle und Dielek- 
trika, endlich durch weitere Erhéhung der Schichtzahl, 
kénnen die Filterkurven natiirlich weitgehend abge- 
aindert werden. So kann man die Linien etwa _ ver- 
schairfen, indem man die Ordnung der Elementarfilter 
erhoht. Wie man durch Anderung der Metallschicht- 
dicken die Linienstruktur andern kann, wurde bereits 
in Abschnitt 4 erdrtert. 

Herrn Prof. Dr. Wotrer danke ich fiir die Anregung 
zu der Arbeit und fiir anregende Diskussionen. 
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Distribution des unités sensorielles dans l’extrafovea * 


J. Ortva et M. AGUILAR 
Instituto de Optica, Madrid 


Sommaire. — Nous étudions ici la symétrie circulaire de la rétine extrafovéale (jusqu’a 28° de la fovéa) en ce qui concerne Vacuité 
visuelle étroitement liée 4 la distribution des unités sensorielles dans la rétine. 
Les résultats obtenus nous montrent que cette symétrie existe et qu’il y a augmentation de la grandeur des unités sensorielles et 
diminution de UV acuité visuelle 4 mesure que nous nous éloignons de la fovéa. 


ZUSAMMENFASSUNG. — Die Arbeit priift die Sehschdrfe und die damit eng verbundene Verteilung der Wahrnehmungselemente auf 
der Netzhaut hinsichtlich ihrer symmetrischen Verteilung auf dem extrafovealen Teil der Netzhaut bis zu einem Abstand von 28° 
von der Fovea. Die Ergebnisse zeigen, dass eine solche Kreissymmetrie besteht und dass je mehr man sich von der Fovea 
entfernt, die Grésse der Wahrnehmungselemente zunimmt und die Sehschdrfe entsprechend verringert wird. 


SumMARY. — A study is made of the circular symmetry of visual acuity of the extra-foveal retina (up to 28° from the fovea), which is 


closely allied to the distribution of sensory units in the retina. 


The results obtained show this symmetry exists, and that there is an increase in the size of sensory units and a reduction of visual 


acuity as the distance from the fovea increases. 


La distribution des cénes et des batonnets dans 
Pextrafovéa a déja été étudiée objectivement par 
OsTERBERG qui trouve une distribution asymétrique. 
Mais ce chercheur ne parle pas des unités sensorielles 
qui nous intéressent particulicrement dans cette étude. 

De nombreux travaux ont été faits sur lacuité 
visuelle et les seuils dans l’extrafovéa. Mais presque 
tous se limitent au méridien horizontal et quelques-uns 
seulement au méridien vertical ; c’est pourquoi nous 
avons pensé intéressant de réaliser une étude sur une 
large zone de l’extrafovéa. 


Observateur;L.C. LUMINANCE :0,6nit 
ORIENTATION SURLA RETINE 


70 10 30 50 70 


Ala 
pe 


10:5 §30° 50 * 76> 10 Moo mcS 70 
ECART ENTRE LESPOINTS (EN MINUTES) 


Fie. 1. 


DISTANCE DE L’ IMAGE A LA FOVEA (EN DEGRES) 


* Ce papier a été présenté partiellement a l’occasion de la 
ata Internationale d’Optique de Florence (septembre 
54). 


135 


180° 


225 


270° 


ee, Ae 


La technique utilisée fut la suivante. 

L’observateur qui regarde deux points lumineux sur 
fond sombre, situés sur une méme horizontale et se 
rapprochant lun de lautre doit signaler le moment 
ou il ne distingue plus les deux points, mais ot il 
apercoit un trait horizontal, et ensuite, le moment ou 
ce trait lumineux se transforme en un point. Toutes 
ces mesures se répétent avec des extrafovéalités diffé- 
rentes ; les résultats obtenus, pour un observateur 
et une luminosité donnés, sont portés sur les courbes 
de la figure 1. 

Sur chacun de ces diagrammes on porte, en ordon- 
nées les extrafovéalités exprimées en grades et, en 
abscisses, l’écart des 2 points exprimé en minutes. 
L’orientation de la vision est indiquée sur chaque dia- 
gramme. Elle est mesurée a partir de la vision nasale 
a laquelle correspond l’angle 0°. La figure 2 fera aisé- 
ment comprendre l’orientation 4 laquelle correspond 
chacune des courbes. Au centre de la circonférence, 
sont placés les deux points constituant le test des 
expériences, et sur la circonférence se déplace la 
lumiére rouge de fixation. L’observateur regarde avec 
Vceil gauche. 
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Les courbes de droite nous donnent la distance a 
laquelle nous cessons de percevoir le trait pour voir 
seulement un point. Ces courbes sont en relation 
étroite avec la grandeur des unités sensorielles, car 
nous croyons que le fait de voir un seul point peut 
faire supposer que les deux images tombent sur la 
méme unité sensorielle. Pour confirmer qualitative 
ment ce fait nous avons réalisé quelques expériences 
en donnant aux points un écart inférieur a celui qui 
permet de voir un trait et en présentant ces points 
alternativement a l’observateur ; celui-ci a vu un 
point constamment immobile au lieu d’un point 
vibrant, ce qui, selon nous, devrait avoir lieu si les 
images des points tombent sur des unités sensorielles 
contigués. En effet, lorsque nous avons fait la méme 
expérience en donnant aux points un écart pour lequel 
ceux-Ci apparaissent comme un trait l’observateur a 
noté une vibration d’amplitude croissante avec l’écart 
des points. 

Les dimensions réelles de ces unités ne coincident 
pas exactement avec les dimensions trouvées expéri- 
mentalement, et cela pour les raisons suivantes : 

1) Les mouvements de image rétinienne pendant 
la fixation de l’ceil, pour une durée de vision de 1/3 de 
seconde (celle employée dans nos expériences) peuvent 
faire varier les résultats d’environ deux minutes au 
maximum (selon Riaes [2]). 

2) Les images des points observés, quoique faisant 
abstraction de la diffraction, ne sont pas ponctuelles, 
mais au contraire, relativement étendues ainsi que 
Va signalé Y. Le Granp [3]. 

Les dimensions des unités sensorielles que nous 
avons trouvées coincident avec celles obtenues dans 
un travail antérieur [4] employant une technique toute 
différente. 

La forme de ces courbes confirme une fois de plus 
VYaugmentation de la grandeur des unités lorsque 
grandit l’extrafovéalité. 

Quoique apparemment ces résultats ne coincident 
pas avec ceux de CLEMMENSEN (1944) [5], selon les- 
quels, les unités sensorielles sont pratiquement égales 
a partir de 3° de la fovéa, en réalité il n’y a pas con- 
tradiction, bien que les unités sensorielles auxquelles 
nous nous référons puissent étre ceiles que nous avons 
appelées unités sensorielles secondaires dans des tra- 
vaux antérieurs [4] ; dans chaque cas notre interpré- 

tation coincide avec celle de CLEMMENSEN lorsqu’il 
' dit qu’a partir de 3° Ja sommation se vérifie grace aux 
chemins transversaux de la rétine, permettant l’inter- 
connection d’un nombre supérieur de récepteurs. 

Les courbes de droite (fig. 1) nous donnent Ja valeur 
inverse de l’acuité visuelle et leur allure nous indique 
la variation de l’acuité visuelle avec lextrafovéalité. 

Sur la figure 3, les ordonnées indiquent la sépara- 
tion des points en minutes, les abscisses, les différentes 
orientations de la vision et dans la partie supérieure 
on exprime en degrés l’angle d’extrafovéalité. 

La courbe inférieure correspond au passage de la 
vision du point a celle du trait, et la courbe supérieure 
au passage du trait 4 deux points séparés. 
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Observateur .L.C. LUMINANCE : 0,6nit 
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Ces courbes montrent clairement la symétrie circu- 
laire existant autour de la fovéa, au moins jusqu’a 
28°, pour ce qui concerne l’acuité visuelle. Cependant 
en ce qui concerne la symétrie nous avons seulement 
expérimenté jusqu’a 28°, nous avons désiré déter- 
miner pour deux des observateurs et avec une seule 
orientation (vision nasale) la prolongation des courbes 
de la figure 1, jusqu’a 60°. Les résultats correspondants 
sont donnés sur la figure 4, et l’on peut voir que le 
comportement des unités sensorielles et de lacuité 
visuelle est analogue 4 celui qui pourrait étre prévu 
par extrapolation des courbes déja obtenues. 

Dans tout ce qui précede nous nous référons a des 
expériences ou les points lumineux étaient situés sur 
la méme horizontale et par suite lorsque nous parlons 
de la variation de grandeur, nous nous référons a la 
dimension horizontale de l’unité sensorielle. Pour véri- 
fier s'il y avait des variations verticalement, nous 
avons répété quelques-unes des expériences avec les 
points situés alternativement en position horizontale 
et verticale, et nous avons constaté que l’écart des 
points nécessaire dans les deux positions était le 
méme. Ceci nous fait supposer une forme a peu pres 
circulaire de ces unités. Bien que ces résultats sem- 
blent contredire ce que l’on sait de l’astigmatisme 
nocturne (6) qui est en général fonction de l’orienta- 
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tion de la lumiére de fixation, cette contradiction est 
Observateur E.B LUMINANCE 0,6nit seulement apparente, parce que le phénomene d’astig- 
ORIENTATION SUR LA RETINE matisme se manifeste dans des expériences ou inter- 
ie) vient la sommation et ce n’est pas le cas ; dans ces 
expériences ce qui influe est la grandeur des unités 
sensorielles, ainsi que la densité de distribution de ces 
unités. 

En résumé nous pouvons dire que 

1) les unités sensorielles dans lextrafovéa ont une 
forme approximativement circulaire ; 

2) la grandeur de ces unités varie en fonction de 
Vextrafovéalité dans toute la zone étudiée (jusqu’a 
; wee 60°) ; 

10 30 = 50 10 2 a 3) cette variation de grandeur garde une symétrie 
ECART DES POINTS EN MINUTES circulaire autour de la fovéa ; 


DISTANCE DE LU IMAGE ALA FOVEA 
N 
@ 
aa | 


4) la densité des unités sensorielles diminue lors- 
Observateur.M.J. LUMINANCE.O,18nit que l’extrafovéalité augmente, et cette densité se 
ORIENTATION SURLA RETINE répartit aussi symétriquement par rapport a la fovéa. 


 Dispositif employé. — La maniére dont s’obtien- 
nent les deux sources est indiquée sur la figure 5. Un 
verre opale de grande homogénéité est éclairé par un 
systeme d’illumination de KOHLER qui permet d’inter-_ 
caler sur Je parcours parallele des rayons, les filtres 
nécessaires pour obtenir Vintensité désirée. On place 
devant le verre opale, une plaque munie d’une fente 
horizontale mesurant 1 mm de largeur. Cette plaque 
est pourvue de plusieurs coulisses lui permettant 
le glissement sur une deuxiéme plaque munie de 


DISTANCE DE UIMAGE A LAFOVEA 


40 60 80 100 120 deux fentes de méme largeur que la précédente et en 
ECART DES POINTS EN MINUTES forme de V, de sorte que l’ensemble des deux plaques 
ne laisse passer la lumiere que par les deux points 

Fic. 4, @intersection de la fente horizontale avec les deux 


PLAQUE PORTANT LES RAINURES EN V autres : il s’agit des deux points lumineux dont l’écart 
PLAQUE PERCEE D'UNE Sw se modifie aisément par le déplacement vertical de la 
FENTE HORIZONTALE deuxieme plaque sur la premiere. 


On a calculé l’angle des fentes en V, afin que deux 
millimetres de déplacement vertical correspondent a 
une variation de deux minutes pour langle qui sous- 
tend les deux points lumineux. On lit le déplacement 
sur une regle graduée en millimetres solidaire de la 
plaque qui porte les fentes enV. _ 


Pour obtenir une lumiére de fixation dont la posi- 
VERRE OPAL tion puisse varier commodément, on dispose une série 
SECTEUR TOURNANT de petites lampes ‘convenablement diaphragmeées et 
pourvues de filtres rouges, sur une baguette rigide 

Fie. 5. mais légére pouvant tourner autour d’un axe perpen- 


/ 
FILTRES 
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diculaire représenté par la ligne horizontale joignant 
Poeil de Pobservateur au point: central des deux lu- 
miueres. hs 
Les lampes rouges sont placées de telle sorte que le 

passage de l'une a la suivante représente une variation 
de langle de fixation de 2°. On peut arriver a 30° de 
la fovéa. Un index permet d’orienter la baguette dans 
des directions formant entre elles 45°. 

. Une résistance variable permet d’établir une inten- 
sité lumineuse convenable. 

‘L’observateur est placé a 2,60 m des points lumi- 
neux, som menton repose sur un appareil dont |’orien- 
tation peut varier, afin que sa face soit toujours diri- 
gée vers la lumiére rouge. 


Mode opératoire. — On choisit premiérement une 


intensité des points qui ne soit pas trop forte et ne 
produise pas une sensation d’éblouissement a |’obser- 
vateur en vision extrafovéale et adapté a l’obscurité 
depuis 40' minutes. L’entrainement des observateurs 
fut laborieux pendant les premiers jours et méme les 
premieres semaines, la précision obtenue fut faible. 

On a amélioré la précision en éclairant les points 
lumineux par intermittence : un secteur rotatoire est 
placé immédiatement devant les points qui apparais- 
sent et disparaissent alors simultanément. Le temps 
d@éclairement fut de 1/3 de seconde, avec un inter- 
valle de 1,7 seconde. 

A part les 40 minutes d’adaptation 4 l’obscurité, les 
séances ont duré approximativement une heure. Chaque 
jour, au commencement, on a situé la lumiére rouge 
dans les positions 16° et 4° de la fovéa sur le cété 
temporal, de sorte que ces deux lectures ont rendu 
possible la comparaison entre les différents jours dont 
quelqués-uns furent exclus par suite de différences 
notables. 

L’intensité des points lumineux a été de 6 x 10°’ 
candelas dans la premiére série d’expériences et de 
1,8 x 10-7 candelas dans la seconde. 

On a-obtenu chaque point par la moyenne de quatre 
lectures réalisées en quatre jours. 

Pour réaliser ces expériences, on a varié le plus pos- 
sible l’ordre d’exposition afin d’éviter que la fatigue 
puisse influencer les résultats finaux. 

Nous avons parlé d’une confirmation qualitative de 
ce que lorsque les deux points apparaissent comme 
n’en formant qu’un seul, leurs images se forment sur 
une méme unité sensorielle. 

Nous avons réalisé cette expérience de la facgon 
suivante. Devant les points et trés pres de ceux-ci 
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on a placé un disque giratoire cachant alternativement 
un point puis l’autre, de sorte que l’apparition de l’un 
d’eux correspondait a la disparition de l’autre. L’ob- 
servateur regardait extrafovéalement et auparavant 
on plagait les points suffisamment prés pour que l’on 
en percoive seulement un. Pendant la rotation du 
disque on ne notait aucune vibration du point lumi- 
neux. Au contraire, lorsque la séparation entre les 
points était telle que Vobservateur voyait un trait 


‘lumineux, pendant Ja rotation du disque, on notait 


une vibration du point. L’amplitude de cette vibra- 
tion devenait plus grande si la séparation des points 


- augmentait. 


Nous avons également fait des expériences pour voir 
quelles modifications se seraient introduites dans les 
valeurs des séparations angulaires entre les points, 
si Pon plagait ceux-ci sur une verticale au lieu d’une 
horizontale. 

Pour cela nous avons monté un dispositif ot lon 
obtenait les points pour un procédé analogue mais en 
substituant les deux fentes lumineuses, en V par deux 
séries de petits orifices d’égale grandeur, disposés éga2 
lement en V. Cette plaque rectangulaire pouvait: 
glisser sur une autre qui portait une fente coupant les’ 
deux rameaux du V afin qu’en chaque position appa- 
raissent seulement deux points lumineux. Le systeme 
était tournant et l’on passait ainsi facilement de la 
position horizontale a la position verticale. 

Nous avons opéré en vision extrafovéale, temporale 
et frontale de 8° et 20°. Nous avons vu que ces valeurs 
obtenues pour la séparation des points étaient indé- 
pendantes de la position horizontale ou verticale de 
ceux-ci, aussi bien quand il s’agissait de voir la sépa- 
ration maximum pour que les points se voient comme 
un seul, que lorsqu’on déterminait la séparation maxi- 
mum des points qui permet de les voir confondus en 
un trait lumineux. 
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Visual experiments with a vibrating test object (*) 


Anna Maria Ercotes. ADRIANA FIORENTINI and G. TORALDO DI FRANCIA 
Istituto Nazionale di Ottica-Arcetri-Firenze 


Summary. — In order to study the influence of eye movements on vision, detection and resolution have been measured with a vibrating 
test object. The experimental data obtained by three observers have been plotted both against frequency and amplitude of the vibra- 
tion. It seems possible to conclude that eye movements, even of a very high frequency, have a réle in affecting vision of moving 


objects. 


Sommatre. — Afin d’étudier V influence des mouvements de l’eil sur la vision, on mesure la limite de perception et le pouvoir de 
résolution avec un test objet vibrant. Les données expérimentales obtenues par trois observateurs ont été représentées en fonction 


de la fréquence et de l amplitude de la vibration. 


Il semble permis de conclure que les mouvements des yeux, méme de trés grandes fréquences, jouent un réle dans la vision d’objets 


en mouvement. 


ZUSAMMENFASSUNG. — Um den Einfluss der Augenbewegungen auf das Sehen zu untersuchen, wurden die Erkennbarkeit und die 
Auflésung an schwingenden Testobjekten gemessen. Die Versuchsergebnisse, die man an drei Beobachtern erhalten hat, wurden 
beide aufgetragen gegen die Amplituden und Frequenz der Schwingung. Man darf wohl schliessen, dass die Augenbewegungen 
selbst bei sehr hohen Frequenzen fiir das Sehen bewegter Objekte eine Rolle spielen. 


1. Introduction. — A question at present of great 
interest in the field of physiological optics is the study 
of involuntary movements of the eye. A number of 
researchers have recently succeeded in measuring 
eye movements, by different methods [1], [2], [3], 
[4], [5], [6], [7]. 

There is some evidence that eye movements have an 
influence on vision. Many interesting researches on 
this subject can be mentioned [8], [9], [10] ; they 
prove that the visibility of very fine objects and the 
resolving power of the eye are greatly reduced when the 
effect of eye movements is eliminated by maintaining 
the image of the test object fixed on the retina. 
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Fic. 1. — Experimental arrangement. 


Some information regarding the influence of both 
voluntary and involuntary eye movements on vision 
can be obtained by another method. Visual acuity 
can be measured with a moving test object. The test 


(*) The research reported in this paper has been made 
possible through the support and sponsorship of the Air 
Research and Development Command, United States Air Force, 
vas its European Office under contract No. AF 61 (514)- 

Some preliminary informations about this research have 
been given at the Florence Meeting on « Problems in Contem- 
porary Optics », 10-15 september 1954. 
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object may be a grating, with black and white parallel 
lines, which is made to oscillate in the direction per- 
pendicular to the lines. Frequency and amplitude of 
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the oscillation is varied in the same range as those of 
natural eye movements. 

The oscillation of the test object produces a displa- 
cement of the image on the retina which combines 
with the displacement due to the eye movements. Of 
course, the type of movements which are easiest to 
realize are regular and periodical, whereas eye move- 
ments are more or less irregular. It is impossible to 
predict the resultant movement. It is possible howe- 
ver, to ascertain by an experimental research whether 
or not visibility of the test object is affected by the 
movement; and possibly what values of amplitude and 
frequency give rise to the greatest effect. This investi- 
gation is the subject of the present research. 
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2. Experimental method. — In order to measure! 
visual acuity with a vibrating test object, the follow- 
ing apparatus has been used. 

The photographic grating R (fig. 1) with parallel 
and alternately opaque and transparent lines, is illu- 
minated through an opal glass by the lamp L, fed with 
d. c. An image of R is focused by the lens Ob very 
close to the plane mirror S. The observer looks at this 
image, reflected by a second mirror S", through the 
semi-transparent plate P. 

Mirror S turns around an axis which makes a small 
angle with the normal to the mirror. During rotation 
of S, the image of R undergoes a circular motion, with 
the same frequency as S. Frequency of rotation is 


and (b) vertical grating. Observer A. F. 
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measured and held constant during each observation 
by means of a stroboscope. The angular amplitude of 
the movement can be adjusted by varying either the 
inclination of S or the distance between S and the ima- 
ge of R (by means of small displacements of mirror S’). 

The component of the circular motion of the image 
parallel to the lines of R has obviously no influence on 
resolution. We must only be concerned with the oscil- 
ation perpendicular to the grating lines. 

If the frequency of the movement were infinite, 
the vibrating test object would be equivalent to a 
fixed test object, with smaller contrast. The luminance 
curve of the equivalent test object can be worked out, 
for any given amplitude of the movement. The compu- 
tations have been performed for the particular case 
used in our experiment where the opaque stripes of 
the test object are three times wider than the trans- 
parent ones. Their photometric profile has been assu- 
med to be rectangular. 
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Figure 2 shows the luminance distribution of the 
motionless test object and of the equivalent test 
objects corresponding respectively to amplitudes of 
1/4, 1/2, 3/4, 1 and 2 times the grating spacing p. 
The relative luminance is plotted against an abscissa 
on the grating taken perpendicular to the lines. 

If now eye movements are taken into account, one 
may anticipate that visual acuity when measured 
with a vibrating test object would not be the same 
as when measured with a test object at rest having 
the corresponding luminance distribution. Indeed it 
is possible as said above that eye movements may 
combine with the test object movement to alter 
visibility. 

The observations have been performed in the follow- 
ing manner. An operator adjusts the parallel line 
test object in a vertical or horizontal position. Then 
he slowly raises the illumination of the grating from 
zero lo a fixed value. (Approximate luminance of 
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— Curves of detection and resolution as functions of oscillation frequency for (a) horizontal 


and (b) vertical grating. Observer M. J. 


dark and bright lines 10 nit and 420 nit). The purpose 
of this procedure is to avoid the effect of an abrupt 
formation of the image on the retina. The observer, 
after a fifteen minute dark adaptation, pulls mirror S” 
towards himself until he is able to detect the orien- 
tation of the lines. The inverse of the angle subten- 
ding the ruling spacing p, for this position of S” has 
been termed the « detecting power » (DP). 

Next, by pulling S” still closer, the observer can 
determine a new position where he is able to resolve 
the single lines of the grating. The inverse of the cor- 
responding visual angle is the resolving power RP. 

The two measurements of DP and RP have been 
performed, for a given amplitude of the movement, 
with each of the following set of vibration frequencies : 
0, 10, 20, 30, 40, 51, 59, 69, 82, 88, 103, 110, 117, 133, 
144 Yr. p.s. 

The succession of different frequencies followed a 
random order, which was unknown to the observer. 
Each set of measurements consisted of two observa- 
tions for each of the above frequencies. Ten sets of 
observations have been made for each of the five 
values of the movement amplitude, namely 1/4, 1/2, 
3/4, 1 and 2 times the grating spacing. Each set of 
measurements took place in a different day. 


3. Results. — Figures 3, 4, 5 show the curves obtai- 


ned by three observers (A. F., A. M. E., M. J.). The 
inverses of both DP and RP (visual angle subtended 
by the grating spacing, in milliradians), are represen- 
ted in the ordinates. The abscissae represent the 
frequency of the oscillation in r. p. s. Each point is the 
mean of 10 independent observations, whose standard 
deviation is of the order + 0.01 milliradian. 

Each curve corresponds to a constant amplitude. 
Amplitudes a for the horizontal test object are respec- 
tively 1/4, 1/2, 3/4, 1, 7/4 times the spacing p of the 
grating. For the vertical test object a has the values 
1/1/4, 2, 3/4, 1, 2 times p. Observer M. J. has not per- 
formed measurements with an amplitude of 1/4 the 
spacing. 

From these curves it is an easy matter to obtain 
other graphs showing DP and RP as functions of the 
amplitude at constant frequency (fig. 6, 7, 8). As 
it is advisable to exclude from our data the accidental 
variations of visual acuity due to a change of the phy- 
siological state of the observer (see the next paragraph) 
we have represented on the vertical axis the difference 
between 1/DP and 1/RP as measured with the moving 
test object and with the fixed one. 


4. Discussion of the results. — When trying to 
interpret the results one should take into account 
that, due to the duration of the observations, the 
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various curves correspond to measurements perfor- 2) Both decreases are fairly regular when fre- 
med at different times. quency inereases from zero to 20-40 r. p. s. 


These measurements are consequently affected by 
the unavoidable changes of the physiological conditions 
of the observer. The value of RP may vary even during 
the same observation. These facts affect in particular 
the values of DP and RP measured with the test 
object at rest. In some case DP or RP for an amplitude 
equal to 1 spacing are smaller than DP and RP for 
an amplitude equal to 2 spacings. The equivalent 
photometric profile has the same shape in both cases 
(fig. 2), but the contrast is smaller when the amplitude 
is twice the spacing : we should consequently expect 
that, at least for the highest frequencies, DP and RP 
would be smaller. This is not confirmed by the expe- 
riments, but we cannot rule out the possibility that 
our result is rather due to a variation of the physiolo- 
gical conditions of the observer than to a different 
effect of wider movements on vision. As regards in 
particular observer A. F. (fig. 3 )) it must be noticed 
that during observations with an amplitude of 1 spa- 
cing, this observer reported a great difficulty to detect 
and resolve vertical gratings. 

For these reasons any conclusion derived from the 
experimental results must be accepted with some 
caution. However it seems that the following points 
can be established with a high degree of reliability. 

1) Both DP and RP measured with an oscillating 
lest object are smaller than those measured with the 
test object at rest. The decrease of DP is less pronoun- 
ced than that of RP. 


3) In about 40% of the curves the decrease of RP 
and DP seems to be followed by a little increase, when 
the frequency becomes a little greater than 20-40 r. p. 
s. A frequency of about 20-40 r. p. s. seems to cause a 
minimum of visual acuity. 

4) When the frequency is greater than 50 r. p. s. 
both DP and RP oscillate about a value approxima- 
tely constant. But several curves, especially those of 
observer A. M. E., show a general decrease also in the 
range of the highest frequencies. 

5) The curves are more irregular in the range of 
high frequencies. 

6) Curves corresponding to vertical test objects are 
more irregular than those corresponding to horizontal 
test objects. 


Finally the curves of DP and RP as functions of 
amplitude of oscillation (fig. 6, 7, 8) suggest the 
following remarks. 

7) For frequencies greater than 10 r. p. s. visual 
acuity always decreases when the amplitude varies 
from 1/4 to 1 spacing, whereas for greater amplitudes 
its behaviour becomes a little uncertain. 

8) When the frequency is 10 r. p. s. visual acuity 
steadily decreases throughout the whole of 
amplitudes. 


range 


It is evident from remarks 2 and 4 that the law of 
decreases of DP and RP undergoes a change about the 
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frequency of 40-50 r. p.s., which is the range of the 
flicker fusion frequency. This fact suggests that in 
the range of low frequencies both eye movements and 
flicker have a rdle in determining the experimental 
values of DP and RP. But, as said in remark 4, many 
curves show a steady though small increase, even in 
the range of the highest frequencies. [t seems probable 
therefore that eye movements of high frequency may 
have also an influence. 

The facet that the movements of 30-40 r. p. s. 
possibly give a minimum of visual acuity (observa- 
tion 3) agrees with the results of CoERMAN (11). He has 
investigated variations of visual acuity of a subject 
during rapid vibrations of the chair on which he was 
seated. The vibrations of the chair are transmitted to 
the subject’s head, and they affect vision in the range 
of frequencies between 16 and 140 ¢. p. s. COERMAN 
found two critical frequencies, giving a minimum 
visual acuity. These le between 25-40 ¢. p. s. and 60- 
90 ¢. p. Ss. respectively. 


AND G. TORALDO DI FRANCIA [Orr. Acra 
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Lettre 4 l’éditeur 


Au sujet de l’aberration sphérique de 1’cil 


ALEXANDRE 
Muséum National d’Histoire Naturelle. 


La remarque de G. Wesrurimer [1] concernant 
Porigine des coordonnées de mes courbes d’aberration 
sphérique de l’ceil est trés judicieuse. Par contre ses 
déductions me semblent exagérées. 

Le fait que dans mes courbes l’aberration sphérique 
tend vers zéro en méme temps que la distance a l’axe 
achromatique, implique en effet que les rayons infini- 
ment voisins de l’axe le coupent au niveau de la rétine, 
autrement dit que la pointe de la caustique se forme 
juste sur la rétine. G. WestHerMeR fait remarquer a 
juste titre que c’est la une supposition gratuite, étant 
donné que nes mesures se limitent a des distances a 
Paxe achromatique supérieures ou égales a 0,5 mm, 
et que par conséquent on ignore la valeur de l’aberra- 
tion pour des distances a axe inférieures a cette 
valeur. I] remplace mon hypothése par une autre, 
a savoir que la valeur de Vaberration sphérique reste 
pratiquement constante pour des distances a laxe 
inférieures a 0,5 mm. Dans ces conditions, aberration 
sphérique de l’ccil moyen serait sensiblement nulle et 
indépendante de l'état d’accommodation, et mes 
mesures correspondraient simplement a un défaut de 
mise au point variable avec effort d’accommodation, 
et autorisé par les conditions expérimentales (test 
constitué de deux fils verticaux, partie utilisée de la 
pupule réduite a deux petits disques). 

Je ferai remarquer tout d’abord qu’il est dangereux 
de raisonner sur une moyenne, surtout lorsqu’elle n’a 
pas grand sens, comme c’est le cas pour l’aberration 
sphérique de l’coil, tres variable suivant lPcoil étudié, 
et méme suivant le méridien considéré [2]. D’ailleurs, 
s'il existe quelque raison de rapprocher aberration 
géométrique de l’oril de l’'aberration dite sphérique [3], 
ceil n’en est pas moins un systéme non centré, et son 
« aberration sphérique » n’est pratiquement jamais 
symétrique par rapport a l’axe achromatique. Si l’on 
examine séparément les résultats que j’ai obtenus sur 
chaque ceil, et des deux cdtés de Vaxe, on verra qu'il 
est impossible de tout interpréter par un simple défaut 
de mise au point, et que l’hypothése de G. WESTHEI- 
MER (aberration constante pour des distances a axe 
inférieures a 0,5 mmr) conduit souvent a des contra- 
dictions de part et d’autre de l’axe. 

Par ailleurs, trois faits expérimentaux au moins 
semblent indiquer que l’aberration sphérique de cer- 
tains yeux continue a varier pour des distances a 
Vaxe inférieures a 0,5 mm : 


1) En utilisant une toute petite région de la pupille, 
on a trouvé que le punctum remotum et le punctum 
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proximum d’un ceil étudié varient avec la distance 
a Vaxe, en sens contraires l’un de l’autre, et que cette 
variation continue nettement pour des distances a 
VPaxe inférieures 40,5 mm [4]. 

2) Parallélement, lorsqu’on a obturé une partie 
axiale de la pupille de 1 mm de diamétre seulement, 
amplitude d’accommodation de l’oril étudié a dimi- 
nué [5], ce qui indique que cette amplitude d’accom- 
modation est particuhérement grande pour les rayons 
paraxiaux (le passage de l’oeil de la sous-correction a 
la sur-correction, joint a la constriction pupillaire se 
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vraisemblablement a certains yeux d’économiser leur 
effort d’accommodation). 

3) En diminuant le diamétre de la région utilisée 
de la pupille, on a trouvé sur deux yeux que l’hyper- 
métropie continue d’augmenter rapidement pour des 
valeurs de ce diamétre inférieures a 1 mm [6]. 


Je pense donc que, malgré certaines observations 
qui tendraient a faire croire le contraire, aberration 
sphérique de I’ceil varie souvent pour des distances a 
axe achromatique inférieures a 0,5 mm. Toutefois 
je concéde tres volontiers a G. WESTHEIMER que cette 
variation est exagérée sur mes courbes, laberration 
nayant pas obligatoirement une valeur nulle sur 
Yaxe. Cette valeur peut d’ailleurs étre obtenue par 
extrapolation d’aprés mes courbes, avec une précision 
pas trop mauvaise, puisque l’on connait la valeur de 
Vaberration 40,5 mm de part et d’autre de l’axe, ainsi 
que la pente générale de chacune des deux moitiés 
de la courbe. 

La figure représente mes résultats moyens obtenus 
en 1953, en prenant égale a zéro Vaberration correspon- 
dant a une hauteur d’ incidence égale a 1 mm, afin de ne 
pas faire jouer un role particuher a Paxe achroma- 
tique. J’ai tracé en trait discontinu les parties des 
courbes relatives a une distance a l’axe inférieure a 
0,5 mm, obtenues par extrapolation comme il vient 
d’étre dit ci-dessus, et par conséquent un peu incer- 
taines (en fait, la variation de aberration au voisir age 
de l’axe semble étre plus rapide pour certains yeux, 
mais serait sensiblement nulle pour d’autres). 

Le passage de la sous-correction a la sur-correction, 
lorsque Pccil accommode de 3 d, se révéle méme sur ces 
résultats moyens et méme si l’on se limite aux dis- 
tances 4 Paxe supérieures a 0,5 mm. Si l’on examine 
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séparément les résultats relatifs a chaque ceil étudié, 
on observe que le phénomene est net pour 5 yeux, 
imprécis pour 4, et pratiquement inexistant pour un 
seul. Les contradictions entre les résultats obtenus par 
divers auteurs ne sont donc pas trés étonnantes, d’au- 
tant moins que ma méthode fournit les résultats rela- 
tifs A un méridien, tandis que la méthode des pupilles 
annulaires fournit un résultat d’ensemble. La méthode 
de parallaxe que j’ai utilisée n’en reste pas moins la 
plus précise, le critérium d’alignement étant trés sen- 
sible, ’'axe achromatique pouvant étre repéré avec 
précision, et les observateurs jouissant librement de 
leur vision binoculaire, ce qui assure la constance de 
la mise au point. L’imprécision relative aux rayons 
paraxiaux peut étre supprimée en effectuant les mesures 
non plus par rapport a Paxe achromatique, mais par 
rapport a un rayon marginal de hauteur d’incidence 
connue. 

Je remercie G. WESTHEIMER de sa critique tres per- 
tinente, qui m’a permis d’interpréter plus correcte- 
ment mes résultats expérimentaux, et qui probable- 
ment mettra fin aux polémiques soulevées ces der- 
niéres années par l’aberration sphérique de |’ceil. 
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